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Sediment  transport  in  the  overland  flow  on  an  impervious  surface 
has  been  studied  to  gain  understanding  of  this  process  in  shallow 
flow  on  a  nonerodible  boundary.  Such  a  case  occurs  in  stormwater 
flow  over  concrete  and  asphalt  pavements  of  roads,  parking  lots  and 
airports.  The  term  sediment  can  include  soil,  rubber  residue  from 
aircraft  tires,  brake  lining  residue,  paint  flakes,  and  other 
particulate  matter.  Previous  and  current  research  has  been 
reviewed  in  order  to  gain  an  understanding  of  the  process  of 
sediment  transport.  Existing  mathematical  models  have  also  been 
reviewed  in  an  attempt  to  find  one  which  can  be  modified  to  give 
predictions  on  the  transport  of  sediment  in  this  case.  Three 
velocity  distributions  (parabolic,  logarithmic,  and  power)  are 
employed  with  a  varying  sediment  diameter  to  depth  ratio  in  order 
to  account  for  the  effect  of  the  sediment  particle  on  the  velocity 
distribution  profile.  The  Shields  diagram  and  previously 
established  plots  of  the  drag  coefficient  versus  Reynolds  number 
for  a  sphere  are  used  to  approximate  the  case  of  incipient  motion 
of  a  particle  in  steady  uniform  and  nonuniform  shallow  flow.  An 
approximate,  theoretically  based  mathematical  procedure  is  derived 
to  estimate  this  process.  w  Appropriate  simplifications, 
assumptions,  or  transformationsViave  been  made  in  order  to  achieve 
a  solution  to  the  model  equatian.  The  end  result  is  a  set  of 
curves  which  relate  the  parameters  of  this  particular  sediment 
transport  problem  using  a  sediment  transport  efficiency  parameter. 
Actual  laboratory  or  field  data  is  not  available  for  this 
particular  case  of  sediment  transport  and  the  resources  were  not 
available  to  obtain  such  data  during  this  research.  However,  the 
presented  model  provides  expected  results  using  hypothetical  data. 
The  assimilation  of  actual  numerical  values  for  the  efficiency 
parameter  and  a  procedure  to  determine  the  total  available  sediment 
is  left  to  others  to  achieve  by  field  or  laboratory  means. 
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modified  to  give  predictions  on  the  transport  of  sediment  in 
this  case.  Three  velocity  distributions  (parabolic, 
logarithmic,  and  power)  are  employed  with  a  varying  sediment 
diameter  to  depth  ratio  in  order  to  account  for  the  effect 
of  the  sediment  particle  on  the  velocity  distribution 
profile.  The  Shields  diagram  and  previously  established 
plots  of  the  drag  coefficient  versus  Reynolds  number  for  a 
sphere  are  used  to  approximate  the  case  of  incipient  motion 
of  a  particle  in  steady  uniform  and  nonuniform  shallow  flow. 
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is  derived  to  estimate  this  process.  Appropriate 
simplifications,  assumptions,  or  transformations  have  been 
made  in  order  to  achieve  a  solution  to  the  model  equation. 
The  end  result  is  a  set  of  curves  which  relates  the 
parameters  of  this  particular  sediment  transport  problem 
using  a  sediment  transport  efficiency  parameter.  Actual 
laboratory  or  field  data  is  not  available  for  this 
particular  case  of  sediment  transport  and  the  resources  were 
not  available  to  obtain  such  data  during  this  research. 
However,  the  presented  model  provides  expected  results  using 
hypothetical  data.  The  assimilation  of  actual  numerical 
values  for  the  efficiency  parameter  and  a  procedure  to 
determine  the  total  available  sediment  is  left  to  others  to 
achieve  by  field  or  laboratory  means. 
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NOMENCLATURE 

Cross-sectional  area  of  flow,  normal  to  flow 
direction 

Plan  area  of  bed 

Areal  grain  packing,  or  the  ratio  of  the  bed 
area  covered  by  the  grains  to  the  total  bed  area 

Cross-sectional  area  of  sediment  particle 

Area  of  channel  bed  covered  by  sediment  particles 

Area  of  sediment  particle  surface  normal  to  flow 

Surface  area  of  sediment  particle 

Distance  between  particle  center  of 
gravity  and  point  of  support 

Characteristic  Diameter  of  sediment  particle 

Force  of  drag  on  a  submerged  sediment  particle 

Force  of  gravity  on  submerged  sediment  particle 

Force  of  lift  on  a  submerged  sediment  particle 

Force  of  hydrostatic  pressure  on  a  submerged 
sediment  particle 

"Pushing”  force  on  a  submerged  sediment  particle 
caused  by  an  adjacent  sediment  particle 

Force  acting  on  a  submerged  sediment  particle 
in  x-direction 

Submerged  weight  of  the  sediment  particle 
Froude  Number 

Weisbach's  resistance  coefficient 

Coefficient  of  resistance  due  to  rainfall 

Gravitational  constant 

Depth  of  overland  flow 

Lateral  input,  such  as  rainfall 

Surface  roughness  of  channel  bed,  including 


any  "roughness"  caused  by  the  formation  of 
the  sediment  particles,  an  average  over  entire 
pavement  area 

Length  of  the  overland  flow 
Momentum  flux 

Sediment  deposition  parameter 
Manning's  roughness  factor 
Wetted  perimeter 

Length  of  sediment  particle  circumference 
contacting  bed 

Pressure  intensity 

Discharge  of  main  flow 

Volumetric  flow  rate  per  unit  width 

Lateral  flow  rate  per  unit  length 
of  channel 

Volumetric  sediment  transport  rate  per 
unit  width 

Dimensionless  sediment  transport  rate 
Hydraulic  Radius 
Reynolds  Number 

Longitudinal  slope  of  surface  of  flow 

Pavement  slope,  the  one-dimensional  maximum 
gradient 

Friction  slope 

Gradient  of  dissipated  energy 

Size  distribution  of  sediment  particles 

Shape  factor  required  to  relate  sediment 
particle  shape  to  that  of  a  sphere 

Storage  of  sediment  within  control  volume 


Free  surface  width 


Mean  flow  velocity 

Critical  mean  flow  velocity  required  for 
incipient  motion  of  the  sediment  particle 

Mean  velocity  of  lateral  flow 

Mean  velocity  of  sediment  particle 

Free  stream  velocity 

Local  temporal  mean  velocity  of  main  flow 
Local  mean  velocity  of  sediment  particle 
Shear  velocity 
Volume 

Sediment  particle  volume 

Settling  velocity  of  a  sediment  particle 

Longitudinal  coordinate  along  channel 
bottom  direction 

Coordinate  perpendicular  to  x  on  vertical 
plane 

Kinetic  energy  correction  factor 

Momentum  flux  correction  factor 

Momentum  flux  correction  factor  for 
flow  local  to  the  sediment 

Specific  weight  of  water 

Specific  weight  of  sediment  particle 

Boundary  layer  thickness 

Nominal  thickness  of  viscous  sublayer 

Parameter  relating  Ucr  to  drag  force 

Angle  between  bed  and  horizontal  plane 

Pressure  correction  factor 

Dynamic  viscosity  of  water 


Kinematic  viscosity  of  water 

Sediment  transport  efficiency  parameter 

Function 

Density  of  water 

Density  of  sediment  particle 

Surface  tension 

Average  boundary  shear  stress 

Sediment  particle  pickup  parameter 

Angle  of  repose  of  submerged  sediment  particle 

Parameter  universal  to  turbulent  flow 

Sediment  particle  weight  parameter 
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1 .  INTRODUCTION 

The  modeling,  and  hence  control,  of  overland  flow  has 
been  a  widely  studied  engineering  problem  for  many  hundred 
years.  The  study  of  erosion,  deposition,  scour  and  related 
problems  is  almost  as  old.  However,  the  fairly  recent 
environmental  "awakening"  by  the  world's  peoples  is  forcing 
the  engineering  profession  to  study  the  less  visible  and 
more  unknown  aspects  of  overland  flow.  One  such  aspect  is 
sediment  as  a  pollutant  carrier.  Although,  the  problem  of 
sediment  as  a  pollutant  and  the  control  of  sedimentation  by 
stormwater  runoff  has  gained  public  attention  and  regulation 
in  general,  the  stormwater  transport  of  pollutants  for 
special  situations  still  requires  much  attention.  This  is 
especially  true  in  the  case  of  pollutant  transport  by  sheet 
flow  across  parking  lots,  airfields  and  other  paved  areas 
into  groundwater,  streams  and  rivers. 

Numerous  equations  and  methodologies  exist  to  calculate 
and/or  estimate  sediment  transport  caused  by  overland  flow 
on  agricultural  and  construction  sites,  and  sedimentation  in 
rivers,  streams,  pipes  and  channels.  However,  little 
research  has  addressed  the  hydraulics  of  sediment  transport 
by  nonuniform  sheet  flow  over  a  nonerodible  impervious 
boundary,  such  as  a  pavement. 

The  scope  of  this  investigation  is  to  apply  dimensional 
analysis,  physical  theory,  and  previously  established 
knowledge  on  fluid  mechanics  and  sediment  transport  to  the 
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case  of  particle  transport  by  overland  flow  on  a  nonerodible 
impervious  boundary.  The  research  is  limited  to  the  case 
similar  to  shallow  flow  normally  encountered  as  stormwater 
runoff  on  airfield  pavements.  Considering  the  given  case 
and  further  restricted  by  the  lack  of  physical  understanding 
and  field  data,  the  objective  of  this  study  is  zo  determine 
the  form  of  a  particle  transport  equation  which  can  be 
applied  to  the  conditions  of  shallow  flow  over  an  impervious 
surface.  This  transport  is  an  approximate  model  given  the 
simplifications  required  for  development.  Future  research 
is  required  to  further  verify  the  model.  The  model  applies 
to  particulate  matter  only  and  is  not  capable  of  simulating 
transport  of  oils,  fibrous  material,  or  dissolved 
pollutants.  However,  future  models  could  account  not  only 
for  the  transport  of  sediment  but  also  other  pollutants  by 
the  use  of  a  concentration  multiplier.  Those  pollutants  in 
dissolved  form  would  be  accounted  for  by  considering  their 
adsorption  to  the  sediment  particles. 
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2.  REVIEW  OF  RELATED  LITERATURE 

2 . 1  BACKGROUND 

Most  of  the  previous  research  in  sediment  transport  has 
been  confined  to  the  cases  of  agricultural  erosion  and  river 
and  reservoir  sedimentation  and  erosion.  The  case  of 
agricultural  erosion  greatly  differs  from  the  case  studied 
here  mainly  due  to  the  presence  of  an  erodible  boundary  and 
non-plane  flow  (i.e.  rills).  The  flow  depths  and  regime 
encountered  in  river  and  reservoir  erosion  and  sedimentation 
make  previous  research  in  that  area  of  limited  value  for 
sediment  transport  by  shallow  flow  on  a  plane,  nonerodible 
surface.  However,  much  literature  concerning  the  above  two 
cases  has  been  reviewed  and  applicable  ideas  and  methods 
have  been  applied  in  this  research.  The  literature  review 
was  divided  into  two  areas  to  simplify  discussion: 
hydraulics  of  overland  flow  and  the  mechanics  of  sediment 
transport. 


2.2  HYDRAULICS  OF  OVERLAND  FLOW 

A  literature  review  of  relevant  fluid  mechanics  and 
hydraulics  was  accomplished  in  order  to  establish  a  sound 
foundation  for  this  research.  The  texts  of  Chow  (1959)  and 
Shames  (1982)  were  particularly  helpful  in  this  regard. 

The  literature  review  on  hydraulics  was  divided  into 


three  areas:  the  study  of  shallow  open  channel  flow,  the 
description  and  measurement  of  resistance  to  flow,  and  the 
description  and  measurement  of  the  velocity  distribution. 
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2.2.1  Shallow  Open  Channel  Flow 

One-dimensional  shallow  open  channel  flow  can  be 
described  mathematically  by  the  Saint-Venant  equations  which 
were  first  developed  by  Barre  de  Saint-Venant  in  1871.  The 
Saint-Venant  equations  are  known  as  the  shallow  water  wave 
equations  and  can  be  expressed  in  terms  of  discharge  or 
velocity.  These  shallow  water  equations  are  composed  of  the 
continuity  equation  (based  on  the  principle  of  conservation 
of  mass)  and  either  the  momentum  or  energy  equation  (based 
on  Newtons'  second  law  of  motion).  The  equations  apply  for 
both  laminar  and  turbulent  flow.  Researchers  use  this 
combination  of  equations  to  describe  shallow  water  flows, 
including  overland  flows. 

Rouse  and  McNown  (1945)  provide  an  interpretation  of 
the  energy  and  momentum  principles  in  hydraulics. 

Various  researchers  have  provided  discussion  on  the 
differences  between  the  use  of  the  energy  and  momentum 
equations.  Yen  and  Wenzel  (1970)  provide  a  detailed 
discussion  of  both  the  momentum  and  energy  approach  and  the 
correlation  between  the  two.  They  also  provide  a  good 
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review  of  past  research  on  equations  for  steady  spatially 
varied  flow.  The  application  of  the  Saint-Venant  equations 
to  the  modeling  of  the  rainfall-runoff  process  is  discussed 
by  Yen  (1986) . 

Steady  uniform  or  nonuniform  flows  have  been  studied 
and  are  presented  for  the  laminar,  transition,  and  turbulent 
regimes.  Unsteady  uniform  flow  rarely  occurs  in  nature  and 
is  not  discussed. 


2. 2. 1.1  Steady  Uniform  Flow 

The  velocity  of  a  steady  uniform  flow  is  usually 
expressed  by  either  the  Manning,  Chezy,  or  Weisbach 
formulas.  The  Manning  formula  is 

C7=-1 ' 49  R2/3S1/2  (2.1) 

n 

where  n  =  Manning's  coefficient  of  roughness, 

U  =  mean  velocity, 

R  =  hydraulic  radius, 

S  =  slope  of  the  energy  line. 

Given  the  geometry  of  the  flow  and  q  (or  U  and  h) ,  Manning's 
n  can  be  found.  Therefore,  n  is  no  longer  an  independent 
variable  but  could  be  used  to  replace  either  U,  q,  or  h. 
However,  this  is  impractical  to  do  so  and  since  n  varies 
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considerably  with  depth,  its  use  is  not  proper  in  shallow 
flow  conditions.  Therefore,  the  use  of  n  has  no  advantage 
for  the  flow  conditions  studied  in  this  thesis. 

The  Chezy  formula  is 

U=CjRS  (2.2) 

where  C  =  Chezy' s  factor  of  flow  resistance. 

The  Weisbach  formula  is 


U= 


N 


jg/Rs 


(2.3) 


where  g  =  gravitational  constant, 

and  f  =  Weisbach  resistance  coefficient. 

The  state  of  steady  flow  can  range  from  laminar  to 
turbulent.  Normally,  the  relationship  between  the  Weisbach 
f  and  the  Reynolds  number  (Re=UR/v)  is  used  to  express  the 
state  of  flow  resistance  for  uniform  and  nearly  uniform 
flows  (Chow,  1959) .  The  Weisbach  equation  can  be  written  as 

f=- gRS  ( 2 . 4 ) 

U2 


In  the  laminar  flow  range  with  Re  <  500, 
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f= 


24 

Re 


(2.5) 


The  Blasius  equation  expresses  the  f-Re  relationship  for 
smooth  pipes  as  (Chow,  1959) 


f= 


0.223 

i?e0-25 


for  150<.Re<25 , 000 


(2.6) 


The  von  Karman-Prandtl  expression  for  higher  Reynolds 
numbers  for  a  smooth  surface  is 


— =2  log {Re/I)  +0 . 4 

s/1 


(2.7) 


And  for  a  hydraulically  rough  boundary 


—  =2  log  ( -^ )  +2 . 16  for  >  50 

if  i?/if 


(2.8) 


where  k  =  the  Nikuradze  equivalent  roughness  size. 

These  equations  are  often  used  for  open  channel  flow.  From 
Eqs.  2.1  to  2.3,  Weisbach's  f,  Chezy's  C  and  Manning's  n  can 
be  related  by 


_8  C 
f~  sfg 


1  •  49  j^i/6 

Jgn 


(2.9) 


2. 2. 1.2  Steady  Nonuniform  Flow 


The  equation  for  steady  nonuniform  flow  with  lateral 
discharge  or  spatially  varied  flow  is  often  solved  for  the 
slope  of  the  flow  surface  (Yen  and  Wenzel,  1970) .  The 
momentum  equation  is 


dh 

dx 


SQ  Sf+ 


Q 

gA 


( Ua  cos4>-2|iJ7)  -  — 

eos6-tg  g 


dp 

dx 


(2.10) 


and  the  energy  equation  is 


dh 

dx 


sa-ss-SL  ( J^ST.  A  cose)  d« 

0  e  gA  gA  2  g  pL _ 2  g  dx 


cos0- 


aU2 

9D 


(2 . 11) 


where  h  =  depth  of  flow, 

A  =  flow  cross-sectional  area, 

D  =  hydraulic  mean  depth  =  A/ (water  surface 
depth) , 

x  =  longitudinal  direction, 

SQ  =  bed  slope, 

Sf  =  friction  slope, 

Se  =  energy  loss  slope, 
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6  =  angle  of  the  bed  with  horizontal  datum, 

P  =  momentum  correction  coefficient, 

a  =  energy  correction  coefficient, 

q  =  lateral  flow  rate  per  unit  length  of 
channel,  negative  for  outflow, 

Uq  =  lateral  flow  velocity  when  joining 
the  main  flow  with  an  angle  <p. 

If  there  is  no  lateral  flow  and  the  variation  of  p  or  a  with 
x  is  small,  Eqs.  2.10  and  2.11  can  be  simplified, 
respectively,  as 


dh 

dx 


SQ  Sf 


cose-!^ 

gD 


(2.12) 


and 


dh 

dx 


S0Se 


COS0- 


g  U2 
gD 


(2.13) 


2.2.2  Resistance  to  Flow 

Resistance  to  flow  is  characterized  by  three  types; 
skin  friction,  form  resistance,  and  wave  resistance.  Skin 
friction  is  usually  the  dominate  resistance  to  flow  in  open 
channels  over  nonerodible  boundaries.  This  skin  friction 
results  from  the  flow  of  fluid  over  the  channel  bed  and 
around  sediment  particles.  In  the  case  of  an  alluvial 
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river,  form  drag  may  also  become  important  in  the  presence 
of  bed  forms  such  as  ripples,  dunes  or  anti-dunes.  Wave 
drag  may  become  important  in  high  Froude  number  flow,  in 
nonlinear  channel  alignment,  and  in  very  shallow  flow  around 
a  large  particle  where  significant  wave  height  is  formed. 

The  open-channel  flow  equations  contain  various 
parameters  to  describe  the  resistance  to  flow.  Rouse  (1965) 
provides  a  form  for  the  rate  of  the  cross-sectional 
resisting  force  and  states  that  the  resistance  coefficient 
in  Weisbach  form,  f,  is  a  function  of  the  Reynolds  number, 
the  relative  roughness,  the  cross-sectional  shape,  the  non¬ 
uniformity  of  the  channel  in  both  profile  and  plan,  the 
Froude  number,  and  the  degree  of  unsteadiness.  The  first 
three  variables  represent  the  effect  of  surface  resistance, 
the  next  two  represent  the  effect  of  form  and  wave 
resistance,  and  the  last  represent  the  effect  due  to  local 
acceleration. 

The  resistance  coefficients  for  a  steady,  spatially 
varied  flow  are  discussed  and  differentiated  in  a  helpful 
paper  by  Yen,  Wenzel  and  Yoon  (1972) ,  The  friction 
resistance  term  found  in  the  shallow  water  equation  is  the 
friction  slope,  Sf.  This  term  is  found  in  the  equation 
derived  by  using  the  momentum  approach,  whereas  the  gradient 
of  dissipated  energy,  Se,  is  used  in  the  energy  derived 
equation.  The  resistance  coefficient,  f,  is  found  as  a 
function  of  either  the  friction  slope,  Sf,  or  the  gradient 
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of  dissipated  energy,  Se.  For  gradually  varied  flow.  Yen, 
Wenzel  and  Yoon  (1972)  describe  the  friction  coefficient  as 
a  function  of  the  Reynolds  number,  Froude  number,  the 
relative  roughness,  the  cross-sectional  shape,  the  non¬ 
uniformity,  and  the  relative  mass,  momentum,  and  energy 
inputs  of  the  lateral  flow. 

Because  of  the  difficulties  involved,  the  friction 
coefficient  is  often  computed  from  the  Weisbach,  Manning,  or 
Chezy  equations  which  were  developed  for  steady  uniform 
flow.  Experiments  on  the  variation  of  f  due  to  impediments 
to  the  flow  (such  as  grasses  or  clumps  of  soil)  have  been 
studied  by  many.  Examples  include  Chen  (1976)  and  Sadeghian 
(1990).  Podmore  and  Huggins  (1980)  used  laboratory 
experiments  to  develop  a  relationship  between  the  physical 
roughness  of  the  flow  area  to  hydraulic  roughness 
coefficients. 

The  value  of  f  is  also  greatly  affected  by  raindrop 
impact.  Woo  and  Brater  (1962)  gave  the  result  of  an 
experimental  study.  Yu  and  McNown  (1964)  used  the  gradually 
varied  flow  equation  and  data  from  the  Army  Corps  of 
Engineers  Los  Angelos  District  to  estimate  the  effect  of 
rainfall  on  the  Weisbach  resistance  coefficient.  Yen  et  al. 
(1972)  and  Shen  and  Li  (1973)  further  expanded  on  the 
raindrop  effects. 

Engman  (1986)  gives  values  for  many  roughness 
coefficients  and  discusses  how  the  coefficients  include  the 


12 


effect  of  raindrop  impact,  channelization  of  flow,  obstacles 
in  the  flow  path,  and  other  factors.  Rouse,  Koloseus,  and 
Davidan  (1963)  describe  f  as  a  function  of  Reynolds  Number, 
channel  shape,  and  the  ratio  of  bed  roughness  to  flow  depth. 


2.2.3  Velocity  Distributions 

A  wide  open  channel  is  one  in  which  the  width  is 
greater  than  10  times  the  depth  of  flow  (Chow,  1959) .  In 
such  cases,  the  velocity  distribution  in  the  center  of  the 
channel  is  not  affected  by  the  sides  of  the  channel. 
Therefore,  the  velocity  distribution  can  be  considered  as 
two-dimensional.  Such  will  be  considered  the  case  in  this 
work.  As  described  by  Chow  (1959) ,  the  velocity 
distribution  is  generally  nonuniform  over  the  flow  cross 
section.  Therefore,  in  the  one-dimensional  formulation  a 
coefficient  must  be  used  to  correct  for  the  velocity  head 
which  is  usually  greater  than  that  computed  by  U2/2g.  When 
the  energy  principle  is  used,  the  coefficient  is  a,  also 
referred  to  as  the  Coriolis  coefficient.  The  velocity  head 
is  then  computed  as  aU2/2g  and  usually  a  ranges  from  1.03  to 
1.4. 

When  the  momentum  principle  is  used,  a  correction 
coefficient  /3  is  used  which  is  also  referred  to  as  the 
Boussinesq  coefficient.  The  value  of  0  commonly  ranges  from 
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1.01  to  1.12.  For  a  two-dimensional  flow, 

h 

Uq$=fu2  dy  (2.14) 

o 


where  u  =  the  local  velocity  component  in  the 

x-direction, 

and  q  =  h  U. 

For  a  three-dimensional  flow, 

U2Ap=Ju2  dA  (2.15) 

A 

The  momentum  coefficient  is  also  defined  and  discussed  in  a 
more  rigorous  manner  by  Yen  (1973) . 

2 . 2 . 3 . 1  Laminar  Flow 


The  velocity  distribution  in  a  steady  uniform  laminar 
flow  is  usually  described  by  a  parabolic  equation 


jL  =  (lz-yl) 
h  h2 


(2.16) 


where  u  =  point  velocity  of  the  flow, 

uh  =  velocity  at  the  water  surface, 
y  =  distance  normal  to  bed, 
h  =  flow  depth. 

The  mean  velocity  is  found  as 
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(2.17) 


2. 2. 3. 2  Turbulent  Flow 

Figure  2.1  describes  the  velocity  distribution  of  a 
typical  turbulent  flow.  As  described  by  Robeson  and  Crowe 
(1985) ,  the  velocity  in  the  viscous  sublayer  can  be 
described  by 


u.  v 


where  u*  =  shear  velocity, 

v  =  kinematic  viscosity  of  the  fluid. 

The  velocity  in  the  turbulent  zone  is  often  described 
by  the  dimensionless,  logarithmic,  universal  velocity- 
distribution  law  (Schlichting,  1968) : 

—  =— (ln-^Hl+JO  (2.19) 

u.  X  v 


where  X  =  a  parameter  universal  to  turbulent  flow  and 

is  independent  of  the  nature  of  the  wall. 

Experimentally,  x=0*4  (Schlichting,  1968).  The  value  of  K 

depends  upon  the  nature  of  the  wall.  The  common  form  of  the 

logarithmic  velocity  distribution  in  turbulent  flow  over  a 
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smooth  boundary  is 

—  =5 .75log  (-H^)  +5.5 
u.  v 


(2.20) 


A  similar  equation  is  used  for  turbulent  flow  over  a  rough 
boundary 


—  =5.75log(^) +8.5  (2.21) 

u.  k 


where  k  =  the  roughness  of  the  channel  bed. 

The  equation  for  a  rough  boundary  is  valid  for  ku./v  >  70. 
Flow  is  said  to  be  hydraulically  rough  when  projections  of 
the  channel  bed  are  fully  exposed  to  the  turbulent  flow. 

This  requires  the  viscous  sublayer  to  be  smaller  than  the 
height  of  the  projections. 

The  region  30  <  yu*/v  <  500  is  the  upper  portion  of 
the  law  of  the  wall  where  the  fluid  viscosity  is  most 
important  in  the  development  of  the  velocity  profile 
(Robeson  and  Crowe,  1985) .  In  the  buffer  zone  where  5  < 
yu./v  <  30  the  velocity  can  neither  be  described  by  the 
viscous  sublayer  distribution  nor  the  logarithmic 
distribution.  It  is  often  assumed  the  velocity  distribution 
of  the  viscous  sublayer  extends  to  the  point  of  nominal 
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thickness  of  the  viscous  sublayer,  <Sn'=ll.  84u/u*  (Robeson 
and  Crowe,  1985) . 

Alternately,  a  power  law  has  been  used  to  represent  the 
velocity  distribution, 

_i 

-  =  (-£)  m  (2.22) 

uh  h 


It  has  been  suggested  that  m  =  6  corresponds  to  Manning's 
formula,  m  =  7  corresponds  to  Blasius'  and  m  =  12  to  the 
formula  of  Hazen-Williams . 


2.3  MECHANICS  OF  SEDIMENT  TRANSPORT 


Sediment  transport,  long  a  topic  of  research  in  regard 
to  erosion,  siltation,  and  river  sediment  yield,  is  becoming 
a  current  topic  of  research  in  regard  to  shallow  water  flow. 
However,  existing  literature  is  not  directly  related  to 
shallow  flow  although  several  approaches  and  fundamental 
findings  are  applicable. 

One  of  the  most  well-known  researchers  in  the  area  of 
sediment  transport  is  A.  Shields  (1936) .  His  plot  (Figure 
2.2)  of  the  dimensionless  shear  stress  versus  the  boundary 
shear  Reynolds  number  is  used  to  predict  the  initiation  of 
sediment  motion.  However,  Shields'  (1936)  work  is  limited 
because  it  does  not  predict  deposition  since  the  critical 
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(from  Graf,  1984) 
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velocity  required  to  pick-up  a  particle  is  different  from 
that  required  for  deposition.  Also,  the  Shields  curve  was 
developed  for  steady  uniform  turbulent  flow  with  a  small 
sediment  particle  size  to  depth  ratio  (ds/h)  .  Its 
applicability  to  nonuniform  and  unsteady  flows  with  a  large 
ds/h  should  be  further  investigated. 

MacDougal  (1933)  used  flume  experiments  at  the 
Massachusetts  Institute  of  Technology  to  derive  a  sediment 
transport  equation  for  bedload  which  varied  with  bed  slope, 
flow  discharge,  and  critical  flow  discharge  at  which 
transport  began.  The  variable  of  mean  bed  shear  stress 
appeared  to  be  the  best  to  describe  sediment  transport  for 
several  flume  experiments  by  Stein  (1965) .  However,  the 
depths  of  flow  in  these  experiments  were  all  greater  than 
0.4  ft. 

Rottner  (1959)  used  dimensional  analysis  to  develop  a 
sediment  transport  model  for  bedload.  The  particle  size  to 
depth  ratio,  ds/h,  was  recognized  as  an  important  parameter. 
Flow  depth  and  velocity  were  also  used  to  compute  the 
discharge;  however,  slope  was  not  included. 

Bagnold  (1966)  approached  the  sediment  transport 
problem  from  general  physics.  Energy  balance  and  mechanical 
equilibrium  considerations  were  used  to  develop  the  model 
which  related  a  sediment  transport  rate  to  a  stream  power 
function.  Yang  (1973)  stated  that  unit  stream  power  is  the 
dominant  factor  to  determine  the  initiation  and  total 
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sediment  discharge  for  an  alluvial  channel  under  equilibrium 
conditions. 

Ackers  and  White  (1973)  used  a  combination  of 
dimensional  analysis  and  physics  to  develop  a  sediment 
transport  equation  for  steady  uniform  flow.  Their 
dimensionless  sediment  transport  parameter  is  described  as  a 
function  of  a  dimensionless  mobility  number  and  a 
dimensionless  grain  diameter.  They  found  that  the  transport 
of  fine  material  was  most  related  to  shear  velocity,  whereas 
the  transport  of  coarse  material  is  most  related  to  mean 
velocity  of  flow.  The  authors  stated  that  their  work  does 
not  apply  for  cases  where  the  Froude  number  is  less  than 
0.8. 

Previous  sediment  transport  laboratory  data  was 
compiled  and  reviewed  by  Cooper  (1970),  Cooper  et  al. 
(1972),  and  Brownlie  (1985).  Since  most  sediment  transport 
formulas  are  empirical  or  semi-empirical  and  developed  from 
laboratory  flume  data,  the  authors  were  concerned  with  the 
dangers  of  engineers  extrapolating  the  equations  to  field 
use.  Most  of  the  laboratory  data  concerned  a  flow  depth  in 
the  range  of  h/ds50  from  10  to  10,000.  The  term  ds50  refers 
to  sediment  particle  diameter  of  which  50%  of  the  particles 
are  finer.  Only  1%  of  the  laboratory  data  was  collected  for 
flows  with  h/ds50  <  7.  Also,  many  previous  data  was 
collected  with  flows  having  a  large  Froude  number  whereas 
the  case  of  shallow,  slow  flow  might  have  a  Froude  number  of 
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0.3.  All  of  the  data  was  collected  for  the  alluvial  channel 
case  where  the  bed  is  erodible  and  the  supply  of  sediment  is 
infinite,  the  Reynolds  number  is  high,  and  the  geometry  is 
not  similar.  Therefore,  the  alluvial  channel  data  is  not 
used  to  verify  the  model  presented  in  this  work  even  though 
some  of  the  data  has  a  ds/h  ratio  in  the  same  range. 

Vanoni  (1984)  provides  a  thorough  review  of  the  history 
of  sediment  transport  study.  Texts  that  are  helpful  in  a 
study  of  sediment  transport  include  Graf  (1984),  Simons  and 
§entiirk  (1976),  and  the  ASCE  Sedimentation  Engineering 
Manual,  edited  by  Vanoni  (1977). 

An  early  study  of  the  effect  of  the  shape  of  a  particle 
on  its  fall  velocity  was  conducted  by  Albertson  (1953). 
Gravel  particles  were  studied  and  the  author  suggested  that 
the  nominal  diameter  of  the  sediment  particle  should  be  used 
when  calculating  Reynolds  number  and  drag  coefficient.  The 
nominal  diameter  is  the  diameter  of  a  sphere  having  the  same 
volume.  This  work  has  been  cited  by  many  researchers 
studying  sediment  transport  in  flows  where  the  depth  is  much 
greater  than  the  particle  size  and  the  particle  terminal 
velocity  is  an  important  parameter. 


2.4  APPLICABILITY  OF  REVIEWED  LITERATURE  TO  THE  STUDIED 
CASE 


Applicable  data  was  not  found  to  calibrate  and  verify  the 
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model  presented  later  in  the  thesis.  Even  data  on  shallow 
flow  with  no  sediment  is  scarce.  The  Shields  diagram  is  the 
only  experimental  data  that  can  be  used  as  an  approximation. 
The  two  most  important  aspects  of  sediment  transport  by 
shallow  flow  are  lower  Reynolds  number  and  large  ds/h  ratio. 
These  two  factors  distinguish  the  case  here  from  that 
generally  studied  previously.  Therefore,  sediment  transport 
relationships  based  on  flows  with  small  particle  size  to 
flow  depth  ratio,  such  as  encountered  in  rivers  and  most 
flume  studies,  are  most  likely  not  applicable  to  this  study. 
Also,  theories  and  ideas  associated  with  flow  of  high 
Reynolds  numbers  are  not  applicable.  This  includes  the 
region  of  constant  Weisbach  resistance  coefficient  at  high 
Re  values. 

Previous  studies  of  sediment  movement  in  turbulent 
flows  have  indicated  that  particles  are  most  often  put  into 
motion  by  the  "gusting"  effect  of  the  turbulent  flow.  This 
is  known  as  the  "bursting  phenomenon."  Since  the  shallow 
flow  is  characterized  by  lower  Reynolds  number,  the  flow  is 
usually  laminar  or  in  the  transition  range.  Therefore,  the 
flow  is  unlikely  to  be  "gusting"  and  the  particles  are  not 
put  into  motion  by  the  "bursting  phenomenon".  Instead,  the 
wake  produced  behind  the  particle  by  the  flow  is  the  primary 
inducement  to  particle  movement.  Since  the  particle  size  to 
flow  depth  ratio  is  large,  the  particles  will  most  likely 
move  by  translation  or  rotation  rather  than  by  lifting  into 


the  flow 
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3.  DIMENSIONAL  ANALYSIS 

The  general  case  of  sediment  transport  by  overland  flow 
is  a  complex  process  which  has  not  been  satisfactorily 
addressed  by  the  engineering  profession.  Rigorous 
analytical  solutions  are  not  available.  Sediment  transport 
by  steady  uniform  flow  has  been  studied  by  several 
researchers  as  indicated  in  the  previous  section.  Because 
of  the  complexity  and  nature  of  the  problem,  all  have  made 
many  and  varied  simplifications  to  achieve  the  final  form  of 
the  model.  However,  many  great  differences  exist  between 
the  previous  work  and  the  focus  of  this  thesis.  The  major 
difference  is  the  lack  of  an  erodible  bed  in  the  case  of 
flow  on  an  impervious  surface  such  as  airfield  pavement. 
Also,  higher  Reynolds  numbers  will  not  be  encountered  in 
shallow  flow  and  the  ds/h  ratio  will  be  much  greater. 
Therefore,  two  approaches  have  been  used  here  in  order  to 
study  the  problem  of  sediment  transport  by  overland  flow  on 
an  impervious  surface.  The  first  approach  is  the  use  of 
dimensional  analysis  and  the  second  is  an  approach  from 
general  physics. 

Figure  3.1  can  be  used  to  represent  the  sediment 
transport  process  for  the  case  studied  here.  A  two- 
dimensional  case  is  considered  and  any  variances  in  the 
cross-section  are  ignored. 
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3.1  DETERMINATION  AND  DEFINITION  OF  SYSTEM  VARIABLES 


The  focus  of  this  study  is  to  arrive  at  an  expression 
for  qs,  the  volumetric  sediment  transport  rate  (per  unit 
width)  for  a  two-dimensional  shallow  water  flow.  Therefore, 
for  the  general  case  of  initially  "clear"  steady  flow  with  a 
very  wide  channel  and  uniform  density,  cohesionless  sediment 
particles  with  no  lateral  input  of  flow  and  no  upstream  or 
lateral  input  of  sediment,  this  rate  can  be  expressed  as: 

qs  =  II1  (m,  P  /  Y  /  °>  Ps»  ds/ SF, SD, Apf  To»  so'  L0' &)  (3-1) 

where  II  represents  a  function.  The  above  variables  have 
been  grouped  into  the  three  categories  of  Fluid,  Sediment, 
and  Flow. 

Fluid  variables  are; 

M  =  dynamic  viscosity, 
p  =  density, 

Y  =  specific  weight, 
a  =  surface  tension. 

Sediment  variables  are; 

qs  =  volumetric  transport  rate  per  unit  width 
(the  dependent  variable  sought) , 

ps  =  sediment  particle  density, 

ds  =  representative  sediment  diameter, 

SF  -  shape  factor, 

SD  =  size  distribution. 
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Ap  =  available  sediment  on  the  bed  expressed  as 

areal  grain  packing  or  the  ratio  of  the  area 
of  the  particles  to  the  total  area  of  the 
channel  bed  in  the  flow  segment  (Note:  the 
number  of  grains  is  a  function  of  A^/ds2)  . 


Steady  uniform  flow  variables  are; 

S0  =  channel  bed  slope  in  the  direction  of  flow, 
h  =  depth  of  flow, 

U  =  cross-sectional  mean  velocity  of  flow, 

)3  =  momentum  flux  correction  factor  due  i  ■* 
a  nonuniform  flow  distribution  over  che 
cross  section, 

L0  =  length  of  flow  segment, 

tq  =  average  boundary  shear  stress 


3.2  SINGLE  PARTICLES  IN  STEADY  UNIFORM  FLOW 


The  sediment  transport  function  of  Eq.  3.1  is  arranged 
into  a  nondimensional  form  with  the  use  of  dimensional 
analysis  and  the  dimensions  of  Length  (L) ,  Mass  (M) ,  and 
Time  (T)  of  the  variables  as  shown  below: 
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Note  that  for  the  case  of  steady  uniform  flow  rQ  = 
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yhS0.  Since  all  other  parameters  are  considered,  rQ  is 
redundant  and  can  be  ignored.  Selecting  p,  h,  and  U  as 
repeating  variables,  Eq.  3.1  can  be  expressed  in  a 
nondimensional  form  as: 


Uhp 


yh 


a 


o  d 

—  ,  —  ,  Sf ,  SDt 
p  h 


.  ,  Ap,  sc,p) 


(3.2 


where  ph\J/n  =  Reynolds  Number  (Re)  , 

yu2p/ (hy)  =  U//gH  =  Froude  Number  (Fr) , 
g  =  gravitational  acceleration, 

U2h/(a/p)  =  Weber  number  (Wb)  . 

Replacing  ps/p  with  (ps/p)-l  and  multiplying  [(ps/p)-l]  by 
the  dimensionless  terms  hUp/ju  and  yh/U2p  gives 

[(Ps/P)-1]  (yh/U2p)/ (M/hUp)  =  [  (ps/p)-l)yh2//xU, 
which  can  be  represented  by  w/U,  where  w  is  the  terminal 
fall  velocity  of  the  particle.  Therefore,  Eq.  3.2  can  be 
expressed  alternately  as: 

=n3  (Re,  Fr ,  Wb,  SF,  SD,  Ap,  ^  ,  S0,  V)  (3.3) 

hU  h  U  v  h 


Many  different  combinations  of  these  parameters  are 
possible.  However,  the  above  are  considered  appropriate  in 
this  study.  Many  of  the  these  combinations  have  been  used 
in  previous  researches.  For  example,  the  shear  Reynolds 


parameter  u.hp/;u  introduced  by  Shields  (1936)  can  be  given 
as  Re/to/U^p  where  u*=>/7071>  is  defined  as  the  bed  shear 
velocity. 
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3.3  DETERMINATION  OF  SIGNIFICANT  DIMENSIONLESS 

PARAMETERS 

Although  all  of  the  above  parameters  apply  to  the  case 
of  transport  of  uniform  density,  cohesionless  sediment 
particles  in  initially  '’clear"  steady  flow  in  a  wide 
overland  surface  or  channel,  a  useable  relationship  cannot 
be  established  with  such  a  large  number  of  parameters.  That 
is,  some  parameters  must  be  combined  or  eliminated  through 
assumptions  or  simplifications  to  reduce  the  total  number  of 
parameters  to  a  manageable  amount.  In  order  to  prioritize 
the  above  nondimensional  parameters,  each  is  summarized 
below: 

Dimensionless  Conditions  When  the  Parameter 

Parameter _  Most  Affects  Sediment  Transport 

Re  The  ratio  of  the  inertia  force 

to  the  viscous  force.  Important 
for  shallow  depth  or  low 
velocities,  i.e.  laminar 
flow  and  transitional  developing 
turbulent  flow  which  is  usually 
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the  case  for  overland  and 
pavement  flows.  Also 
important  to  show  the  effect 
of  temperature  changes  on  the 
flow. 

Fr  Ratio  of  the  inertia  force 

to  the  force  of  gravity. 

Without  surface  wave  distortion, 
it  is  not  important  when  the  flow 
is  subcritical  with  Fr  <  0.5. 

Wb  The  surface  tension  effect  is 

expected  to  be  unimportant  unless 
the  amount  of  water  is  so  small 
that  a  continuous  sheet  flow 
cannot  be  sustained  and  the 
flow  is  broken  into  small  pots, 
which,  if  this  ever  occurs,  it 
will  only  be  at  the  upstream  end 
of  the  overland  surface. 

ds/h  It  reflects  the  particle  size  as 

an  obstacle  to  the  flow.  When 
its  values  is  small,  it  merely 
reflects  the  surface  roughness  of 
the  bed.  When  ds  is  about 
one-quarter  of  depth  or  more, 
it  also  reflects  the  disturbance 


of  the  particle  to  the  flow 
velocity  distribution. 

The  relative  mobility  of  the 
particle  in  water.  Since  the 
flow  depth  is  small,  the 
particle  is  never  able  to 
reach  terminal  fall  velocity,  w. 

In  the  case  of  a  constant  ps/p, 
this  parameter  is  not  significant 
in  shallow  flow  and  w  is  used  only 
as  a  sediment  property  indicator. 
However,  this  parameter  is  useful 
in  the  study  of  nonuniform  flow 
since  it  is  related  to  the  bed 
shear  tq  and  bed  slope  S0. 

The  shape  of  the  particle  affects 
the  drag  force.  Thus,  SF  is  most 
important  when  the  shape  is  far 
from  spherical,  e.g.  a  disk 
normal  to  the  flow.  For  normal 
sand-gravel  type  sediment,  the 
significance  of  this  parameter 
is  only  of  secondary  importance. 
Significant  only  for  greatly 
varying  particle  sizes.  Equal 
to  zero  for  uniform  particle  sizes. 
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L^h  Since  this  term  expresses  the 

length  of  the  flow  channel  where 
sediment  can  be  picked  up,  it 
must  be  considered  for  the 
calculation  of  sediment 
transport. 

Ap  Most  important  when  the  sediment 

density  is  sparse.  For  densely 
packed  particles  such  as  in 
alluvial  channels,  the  sediment 
packing  density  becomes 
statistical . 

SQ  The  channel  bed  slope  reflects 

significance  of  bed  shear,  and 
thus,  the  ability  of  the  flow 
to  move  the  bed  particles. 

P  When  the  velocity  distribution 

is  highly  nonuniform  and  varying 
with  the  flow  distribution. 
Normally,  in  cases  of  a  large 
flow  depth  to  sediment  particle 
size,  the  effect  of  this  parameter 
parameter  is  small.  However,  this 
is  not  the  case  here. 

Using  the  above  evaluations,  for  a  steady  uniform. 
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shallow,  subcritical  flow  on  an  impervious  pavement  of  small 
slope  with  nearly  spherical  constant  density  sediment  on  it, 
the  dimensionless  independent  parameters  in  Eq.  3.3  can  be 
prioritized  into  primary  and  secondary  importance  categories 
as  follows: 

Primary  Parameters  Secondary  Parameters 

Re ,  ds/h ,  Lyh ,  Ap , SQ ,  £  Fr ,  w/U,  SF,  SD,  Wb 

Thus,  Eq.  3.3  can  be  reduced  to 

-^.=11 4  (i?e,  —  ,Apl$,SQI  — )  (3.4) 

If  qs/qs  ^  versus  flow  length,  L,  were  plotted,  the 
maximum  sediment  transport  rate  is  achieved  after  a  finite 
distance  and  does  not  increase  with  additional  downstream 
distance.  Since  qS(nax  would  be  difficult  to  determine,  Lx  is 
defined  as  the  flow  length  corresponding  to  a  sediment 
transport  rate  which  is  x%  of  the  maximum  transport  rate. 

Lx  will  be  referred  to  as  the  threshold  value  of  L,  flow 
length.  Above  a  value  corresponding  to  L^,  the  sediment 
transport  rate  is  almost  independent  of  L0  and  hence  L^h. 

The  actual  value  of  x  could  be  set  as  95,  98,  or  some  other 
meaningful  level. 

A  dimensional  analysis  is  accomplished  in  the  same 
manner  as  that  for  qs  in  order  to  develop  the  function  for 
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Ucr,  the  critical  flow  velocity  required  for  incipient 
motion  of  the  sediment  particle.  Neglecting  the  terms  a, 
SF,  SD,  Ap,  and  Lp,  the  dimensionless  relationship  for 
critical  velocity  is 

J^^(^,^,Rbc,S„2s)  (3.S) 

n  P 


where  Rec  =  (Ucrhp)/*i, 

Ucr/^s  =  (PUcr2)/(Yds). 

For  both  above  cases,  steady  nonuniform  flow  can  be 
represented  by 


(*e,  ,  A  p,  s  is  ,  _I2_ 

hU  5  h  p  0  h  p  u2 


(3.6) 


where  Tp/pU2  is  used  to  represent  the  dimensionless  bed 
shear  and  accounts  for  the  nonuniformity  of  the  flow. 

In  order  to  consider  the  effect  of  rainfall  on  the 
sediment  transport  rate,  other  parameters  must  by  included. 
As  given  by  Yoon  and  Wenzel  (1971),  the  dimensionless 
relationship  for  spatially  varied  flow  is 


uo  _ 


i  rr2 


=nfi  (Re,  Fr, 


p  ih 


.Fr. 


P  idK  e 


,^S0,T 1  ,  X) 


(3.7) 


where  FrR  =  Froude  number  of  the  rainfall, 

i  =  rainfall  intensity, 
c  =  boundary  roughness, 


rj  =  raindrop  pattern  parameter, 

A.  =  raindrop  shape  parameter. 

The  case  of  unsteady  flow  can  be  represented  by  including 
term  to  account  for  the  variation  of  the  flow  depth  with 
time.  The  general  case  would  include  all  of  the  above 
considerations . 
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4 .  THEORETICAL  CONSIDERATIONS 

The  sediment  transport  model  is  developed  by  reviewing 
physical  considerations  and  comparing  the  outcome  of  this 
review  with  the  results  from  dimensional  analysis.  Past 
developments  in  sediment  transport  technology  are  modified 
for  the  case  of  shallow  flow  over  an  impervious  surface. 
Modifications  are  necessary  to  account  for  the  shallow  flow 
depth  to  particle  diameter  ratio  and  also  to  reflect  both  a 
laminar  and  turbulent  velocity  distribution. 

First,  the  particle  mobility  problem  is  addressed  from 
general  physics.  The  first  procedure  is  a  definition  of  all 
forces  acting  on  a  particle  and  the  critical  condition  when 
the  particle  is  about  to  move.  The  second  procedure  is  the 
application  of  Newton's  law  of  motion  to  a  sphere. 


4.1  FORCES  ON  SUBMERGED  SINGLE  PARTICLES 

4.1.1  Steady  Uniform  Flow 

Figure  4.1  depicts  a  particle  at  rest  on  a  plane  bed 
submerged  under  steady  uniform  flow.  Incipient  motion  is 
described  as  the  point  at  which  the  particle  is  just  about 
to  move  if  any  additional  force,  however  small,  is  applied. 
Once  this  additional  force  is  applied,  the  particle  will 
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move  by  ary  one  or  a  combination  of  three  means;  lifting, 
sliding,  or  rotating. 

Forces  acting  on  the  particle  in  Figure  4.1  include 
drag,  lift,  weight,  and  resistance.  The  drag  and  lift 
forces  are  due  to  the  flow  surrounding  the  particle.  The 
weight  force  is  the  submerged  weight  of  the  particle  in  the 
fluid.  Resistance  is  the  friction  force  acting  between  the 
particle  surface  and  the  bed  surface.  The  drag,  lift, 
weight,  and  resistance  forces  can  be  described  by  the 
following  equations: 

Liag:  F 'D- C d~2~  ^  sn  (4.1) 

ifjr 

Lift :  FL=CL-L-Asn  (4.2) 


Weight:  Fv=(ye-y)  Ve=  (ps-p)  gVs  (4.3) 


Resistance:  FI=Cii(Fwy-FL) 


(4.4) 


where  A  =  the  projected  area  of  the  particle  normal  to 
the  direction  of  flow, 

C0  is  the  drag  coefficient, 

CL  is  the  lift  coefficient, 

U  is  a  reference  velocity  which  can  be  related  to 
tne  approaching  flow  velocity,  U,  along  the  x- 
direction  passing  through  the  particle's  center  of 
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gravity,  G,  and  a  coefficient  C^,  i.e., 

Ur  =  C^U2,  Vs  is  the  volume  of  the  particle, 
and  is  a  nominal  friction  coefficient. 

The  coefficient  Cu  is  a  function  of  the  momentum  flux 
correction  factor,  0. 

Lift  occurs  when  the  lifting  force,  F,  ,  is  greater  than 
the  submerged  weight  of  the  particle  in  the  y-direction, 

Fyy.  The  case  of  incipient  motion  by  lifting  is 


Fl-F^O  (4.5  a) 

or 

CL~~Asn=  (P.-P)  <7^C0S©  (4 . 5b) 

Solving  for  the  incipient  motion  critical  mean  velocity, 

Ucr,  provides 

£&=(-^-l)  g-^cosG^  (4.6) 

where  K1  =  2Vs/ASn  2-s  a  function  of  the  geometry  of  the 

particle  and  is  equal  to  8/3  for  a  sphere. 

The  lift  coefficient  CL  is  a  function  of  the  Reynolds  number 

of  the  flow  around  the  particle  (Re  =  Uds/v) .  The 

coefficient  Cu  is  a  function  of  the  momentum  correction 

coefficient  over  the  depth  of  the  particle,  0S, 
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P  s= 


(4.7) 


where  Ur  =  a  reference  velocity  of  the  flow. 

Hence,  for  a  constant  0,  Ucr  can  be  expressed 

UCI=KL(PsfRe)ja;  (4.8) 

where  KL  =  a  function  of  /3S  and  Re. 

For  the  airfield  case,  pavement  slopes  range  from  2%  to  7%. 
Therefore,  e  ranges  from  approximately  1.15  to  4.01  degrees, 
and  cose  ~  1,  sece  -  l,  and  tan©  -  sine  =  S0. 

The  particle  will  slide  along  the  bed  surface  when  the 
resultant  force  applied  by  the  flow  in  the  x-direction  is 
positive.  The  incipient  case  of  translation  is 


Fd+Fw  sin0-Fr=O 


(4.9) 


Substituting  Eqs.  4.1  to  4.4  into  Eq.  4.9  and  solving  for 
the  incipient  motion  critical  velocity, 


r  ( C„  cos0-sin0)  ,  p_  4  ,  .  4 


(4.10) 
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Therefore, 

uCi=kd(  Ps'  Re>  Cy)  (4.11) 


where  =  a  function. 

Even  if  the  particle  does  not  move  by  lifting  or 
sliding,  it  can  still  be  transported  by  rotating  or  rolling 
along  the  bed.  This  would  occur  when  the  forces  exerted  in 
the  x-direction  about  the  support  point  of  the  particle 
create  sufficient  torque  to  overcome  the  resisting  torque 
(Figure  4.2). 

The  incipient  case  is 

FDeD-FLeL+Fu-^  tan6=0  (4.12) 


The  moment  arms  eD  and  eL  are  a  function  of  the  particle 
geometry  and  flow  velocity  around  it,  and  hence,  are  a 
function  of  0S.  Therefore,  by  substituting  Eqs.  4.1  to  4.3 
into  Eq.  4.12,  the  critical  velocity  can  be  expressed  as 

U^K^'.Re)^  (4.13) 

In  each  incipient  case  of  lifting,  sliding,  or  rolling, 
the  critical  condition  is  expressed  in  terms  of  Ucr,  the 
critical  approaching  flow  velocity  in  the  x-direction 
passing  through  the  center  of  gravity  of  the  particle.  The 
main  observation  from  the  above  physical  view  is  that  the 
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critical  velocity,  Ucr,  is  a  function  of  yds,  or  a  function 
of  the  particle  weight  to  the  1/6  power.  This  relationship 
of  critical  velocity  to  weight  is  known  as  the  one-sixth 
power  law  of  A.  Brahms  of  1753. 

Note  that  in  all  cases,  the  term  Asn,  the  particle 
surface  area  in  contact  with  the  flow,  changes  according  to 
the  degree  of  embedment  of  the  particle  as  shown  in  Figure 
4.3.  Therefore,  this  term  alone  can  account  for  all  degrees 
of  embedment . 


4.1.2  Steady  Nonuniform  Flow 

Figure  4.4  depicts  a  single  particle  submerged  in  a 
steady  nonuniform  flow.  Incipient  motion  is  conditioned  in 
the  same  manner  as  stated  for  the  uniform  case.  However, 
the  particle  is  also  subject  to  the  force  due  to  the 
pressure  gradient  of  the  nonuniform  water  surface  profile. 
The  equations  for  drag  and  lift  are  rewritten  in  the 
following  manner 

Drag  -  [  j pr assure  cL4s]  X-Component+  f  J shear  cL4s]  x-component  (4  ■  14a) 

Am  A9 

Lift  -  [  j"  pr  aSSUre  cL4s]  y-component+  t  J"  Shear  dAs]  y.Coirf>onent  (4  •  14b) 
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The  first  integral  represent  the  pressure  drag  and  the 
second  term  represents  the  shear  drag.  The  equation  for 
rotation  is  modified  similarly.  Therefore,  for  flow  of 
small  or  moderate  nonuniformity,  the  critical  flow  velocity 
can  be  written  as 

(4.15) 


As  stated  by  Graf  (1984,  p  37): 

Investigations  have  shown  that  there  does 
exist  an  intricate  relationship  between  the 
inertia  and  viscous  forces,  and  only  under  certain 
conditions  does  the  problem  lend  itself  to 
mathematical  theory. ..." 

Therefore,  most  of  the  complications  of  the  viscous  effects 
are  summarized  in  the  drag  and  lift  coefficients  CD  and  CL, 
which  are  determined  by  experiment  with  the  same  Reynolds 
number  in  the  laboratory  flume  as  that  prevailing  in  the 
field.  However,  most  CD  or  CL  vs.  Re  relationships  for 
spheres  are  for  non-rotating,  smooth  particles  surrounded  in 
a  fluid  with  constant  approaching  velocity. 


4.2  THE  MOMENTUM  APPROACH 
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The  principles  of  continuity  and  momentum  can  be 
applied  to  the  control  volume  indicated  in  Figure  4.5  to 
obtain  a  sediment  transport  equation.  Considering  the  fluid 
and  sediment  as  two  separate  control  volumes  will  produce 
four  equations;  both  a  continuity  and  momentum  equation  for 
the  sediment  and  fluid  volumes.  The  continuity  equation  for 
the  steady  flow  without  lateral  flow  is 


0 

dx  dx 


(4.16) 


The  momentum  equation  for  the  fluid  is  given  as  Eq.  2.12 


dh  _  S0  Sf 


(2.12) 


where  D  =  A/ (water  surface  width)  is  the  hydraulic  mean 

depth . 

The  friction  slope  Sf  can  be  estimated  using  the  Weisbach 
formula  Sf  =  fU2/8gR,  or  the  Manning  formula,  Sf  = 

(n2U2)/(2.22  R4/3)  .  For  unit  width,  the  cross-sectional  area, 
A,  is  equal  to  the  depth,  h.  Also,  in  steady  uniform  flow, 
dh/dx  =  0,  and  hence  the  flow  surface  slope  is  equal  to  the 
i^'ction  slope. 


4.2.1  Sediment  Equations 
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Assuming  a  sparse  concentration  of  sediment  in  the  flow 
on  the  channel  bed,  the  continuity  equation  for  the  sediment 
is  given  as 

Qs2=Qsi  +  ^*x  (4.17) 

where  qs  =  the  sediment  discharge  per  unit  width  of 
channel  in  volume. 

At  the  upstream  end  of  the  overland  surface,  qs1  =  0.  As 
will  be  ;hown  later,  the  sediment  volume  flux  within  the 
control  volume  is 


where  n  =  the  number  of  sediment  particles  on  the  bed  in 

the  control  volume, 

Vs  =  K2ds3  =  the  volume  of  a  single  particle 
l  =  a  sediment  transport  efficiency  parameter. 

The  number  of  sediment  particles  is  a  function  of  the  areal 
packing,  Ap,  and  the  square  of  the  characteristic  diameter 
of  the  particle,  ds.  The  areal  packing  is  defined  as 


a'  nd2sK , 
A=  —  = — 

P  Ab  dl 


(4.19) 


where  As'  =  the  area  of  the  channel  bed  covered  by  the 

sediment  particles, 

Ajj  =  the  area  of  the  channel  bed, 

K1  -  a  function  of  the  particle  shape  and  is  equal 
to  ir/ 4  for  spheres. 

Substituting  Eq.  4.18  into  Eq.  4.17  provides 
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Qe2=Qsl+nZK2dl&x 


(4.20) 


The  momentum  equation  for  the  sediment  volume  can  be 
expressed  as 


5>« 


dmsUs 

dt 


(4.21) 


where  Fxs  =  a  force  acting  on  the  particles  in  the 

x-direction, 

ms  is  the  mass  of  the  sediment  particles, 

U  is  the  velocity  of  the  sediment  particles  in 
the  x-direction. 

Since  the  acceleration  of  the  particles,  dUs/dt  is  generally 
unknown  and  rather  difficult  to  find,  Eq.  4.19  provides  no 
useful  help  in  solving  for  qs. 

Instead,  assuming  for  the  particles  collectively  the 
acceleration  E  (dUs/dt)  =0 ,  the  force-momentum  relationship 
for  the  fluid-sediment  mixture  of  the  control  volume  can  be 
written  as 


-2-  (p^g^  +  pPAtf2)  =VsSQAs+ySQA-TbPs-x0  (P-Ps) 


(4.23) 


where  r0  =  the  bed  shear  over  the  wetted  perimeter  P  not 
covered  by  the  particles, 

rb  is  the  shear  on  the  particles  on  the  boundary 
occupying  a  perimeter  Ps. 

The  bar  above  /9S  and  Us  indicates  averages  over  the 
particles,  and  the  particles  velocity  can  be  related  to  U  by 
Us  =  l^U,  where  the  coefficient  Kj  is  a  function  of  the  flow- 


51 


particle  momentum  correction  coefficient  /3S. 


4.3  DETERMINATION  OF  VELOCITY  DISTRIBUTIONS  FOR  DIFFERENT 
MOMENTUM  CORRECTION  COEFFICIENTS 


The  purpose  of  this  section  is  show  the  relationship 
between  the  point  and  mean  velocities  of  the  flow  over  the 
flow  depth  or  height.  Parabolic,  logarithmic,  and  power 
velocity  distributions  shown  in  Figure  4.6  will  be 
considered. 

Although  the  particle  is  free  to  move  in  a  three- 
dimensional  manner,  only  the  two-dimensional  case  will  be 
considered  here.  Furthermore ,  only  the  movement  in  the 
down-channel  direction  (or  x-direction)  contributes  to  the 
sediment  transport  rate,  qs. 

Three  Boussinesq  coefficients  are  considered  for  the 
study  of  sediment  transport  in  shallow  flow  where  the  flow 
depth  to  sediment  size  ratio  is  small.  One  coefficient,  /3, 
represents  the  correction  factor  for  the  entire  flow  depth, 
h,  and  the  other  two  relate  to  the  depth  at  the  particle 
height,  ds.  These  coefficients  are  defined  by  Eq.  4.7  where 
the  reference  velocity  of  the  flow  is  now  defined  as  the 
mean  velocity  at  the  particle  height  for  0s1 ,  and  the  mean 
velocity  at  the  flow  depth  for  (3  and  0s2.  Therefore, 


Figure  4.6  Definition  sk 
near  Particle 
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(4.24) 


(4.25) 


d. 


$.2 


J  u2  dy 

_o _ 

U2h 


(4.26) 


where  u  =  the  x  component  of  the  local  point  velocity  of 

the  fluid 

and  u1  =  the  average  fluid  velocity  over  the  particle 

depth. 


d. 

u[=-±-fudy  (4.21) 


A  plot  of  (3  and  /?s  versus  ds/h  should  produce  a  useful 
family  of  curves.  For  example,  as  discussed  in  Section 
4.1.1,  the  incipient  motion  critical  flow  velocity  is  a 
function  of  fl  .  Since  the  situation  under  review  concerns  a 

^  S 

wide  channel,  0  and  /?s  versus  ds/h  represents  a  unique 
solution.  This  would  not  be  the  case  for  a  narrow  channel 
where  velocity  distribution  is  modified  near  the  surface  due 
to  the  influence  of  the  walls.  Therefore,  the  three  types 
of  velocity  distributions  mentioned  above  will  be  evaluated. 
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4.3.1  Parabolic  Velocity  Distribution  for  Laminar  Flow 


It  is  assumed  that  laminar  flow  considered  here  can  be 
represented  with  a  parabolic  velocity  distribution  of  the 
form 


u  _  /  2  y_y2) 

uh  h  h2 


(4.28) 


where  h  =  the  depth  of  flow. 

The  above  equation  gives  the  maximum  velocity  at  the  water 
surface,  where  u  =  uh.  Substituting  Eq.  4.28  into  Eqs. 
4.24,  4.25,  and  4.26  and  solving  provides 


P  = 


2  * 


U2h 


/< 


4 y2  _  4y3 
h2  h2 


dy 


(4.29) 


Integrating  Eq.  4.29  and  substituting  Eq.  2.17,  U  =  2uh/3, 
into  it  yields 

|u2P  =  _Lu2,  p  =!  .  200  (4.30) 


Similarly,  0s1  is  derived  by  substituting  Eqs.  4.27  and  4.28 
into  Eq.  4.25  and  /3s2  by  Eq.  4.28  into  Eq.  4.26.  By 
applying  Eq.  2.17  to  the  result  of  integration,  one  has 
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P 


31 


9 

5 


+ 


d~  c  (j  2 

15-10(-f)+4(-5) 
h  3  h 


(4.31) 


^-3<T)  -7(' 


A  7  20 


d  5 

T> 


(4.32) 


Figure  4.7  is  a  graph  of  0s1  and  /3s2  versus  ds/h  for  the 
parabolic  velocity  distribution. 


4.3.2  Logarithmic  Velocity  Distribution  for  Turbulent 
Flow  over  both  Smooth  and  Rough  Boundaries 

The  turbulent  flow  velocity  distribution  has  been 
discussed  in  Section  2. 2. 3. 2.  For  a  logarithmic  velocity 
distribution  over  a  smooth  boundary  shown  in  Figure  2.1, 

-jj-=y.  0  sy.sfi,  (4.33) 

-jj-=5.7  5  log  y,  +  5.5  b'<.y.<Lh  (4.34) 

where  y*  =  yu./v  and  <S.  =  <Su,/v 

S  =  the  nominal  thickness  of  the  viscous  sublayer 
and  is  assumed  equal  to  11.6  for  the  smooth 
boundary. 

For  turbulent  flow  over  a  rough  boundary,  Eq.  4.33  is 
assumed  to  still  be  valid  if  not  obstructed  by  the  roughness 
element  and  approximately  for  <5*  <  y,  <  h 
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-jj-  =5.75  log  y.  -5.75  log  Jc.+  8.5  (4.35) 


where  k*  =  (k  u*/v)  where  k  is  the  roughness  height. 

For  the  rough  boundary  case  <5*  may  be  approximated  as  the 
intersection  of  Eqs.  4.3  3  and  4.3  5  which  implies  <5.  being  a 
function  of  k*.  Accordingly,  it  has  been  found  that  <5. 
ranges  from  2.49  to  11.6  for  which  k.  ranges  from  27.63  to 
3.35. 

Substituting  Eqs.  4.33  and  4.34  into  Eq.  4.24  provides 


P  = 


1 

U2h 


(5.75  logy.+5 . 5)  2  cfyj 


(4.36) 


Integrating  Eq.  4.36 


P=l  + 


(33.06  logh,+17 . 26 ) 2-  517  ~59  +5.95 

_ K _ 

(  U(i)  )  2 
u. 


(4.37) 


where 

— —  =5.75  log  h.+3.003--r-^  (log  6.-3.003-—^)  (4.38) 

u,  a  h.  2 

Similarly,  substituting  Eq.  4.33  and  4.34  into  Eq.  4.25  and 
4.26,  respectively,  and  integrating  provides 
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(32.06  log  dg,  + 17 .26)2---? '59-+5.95 

{  ad,(»)  )2 

U. 


6  =i  +  — S. 
Ps2  J-  h 


d  (32.06  log  dfl.+17.26)2---^--59  +5.95 


U. 


where 


O', 


^^-=5.75  log  d  +3 .003--^  (log  6, +3.003--) 
u,  s  h.  2 


For  5,  =  11.6,  Eqs .  4.37,  4.39,  and  4.40  become 


P-1- 


(32.06  log  A. +17 .26)2--™' l'59- +5. 95 

_ K _ 

A* 

u . 


(32.06  log  d  + 17 .26)  2---1?  •  59  +5. 


95 


PSi  =  l  +  - 


(  — )2 

U. 


P^=1+T 


(32.06  log  d  +17 . 26 )  2-  517 /5i+5.95 
d0  s  d„. 


(^4)2 

u. 


(4.39) 


(4.40) 


(4.41) 


(4.42) 


(4.43) 


(4.44) 
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where 


—£■  =  5.75  log  A.+3.003--38'55  (4.45) 

u,  h. 


-^2=5.75  log  d+ 3.003-  3  8.'  5  5  (4.46) 

U.  <4. 


In  order  to  express  Eqs.  4.39  and  4.40  in  terms  of  ds/h,  the 
following  equations  are  substituted 


(4.47) 


log  c?s,=log 


(-^)  +log  h. 
n 


(4.48) 


log  ds=log ( — )  (^) 


(4.49) 


log  c?s,/c.=log 


(4.50) 


Therefore,  using  Eqs.  4.47  and  4.48,  /3s1  and  (3sZ  can  be 
expressed  as 
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[32.1  (log-^g  +logh.)  2+17 .3] 2-  —  7--6  (-§-)  +5.9  5 
P.,-1* - - - Jj - -  (4.51) 

(Z±)2 

u. 


[32 . 1  (log-^  +logil.)  2 +  17 .26]2-51?-6  +5.95 

.  A  h,  a„ 

P£2=l  + - - - - - 2 -  (4.52) 

(-2)2 

U. 


For  the  rough  boundary  case,  Eqs.  4.33  and  4.35  are 
substituted  into  Eq.  4.24 


P  = 


1 

U2h 


«.  d,. 

[  fy*dy,+  f  (5.75  (log  y.-log  1c.)  +8 .  5)  2c?yJ 
0  ft. 


(4.53) 


Integrating  Eq.  4.53 


6.  63. 

b~Tb'T 

p=l  + - 77 -  (4.54) 

(-1)2 

u. 


where 
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Up  h,  6.  ..  hm  6.  v 

—  =log  -^+1.04-—^  (log  —^+1.04--^) 
u.  k ,  h,  k.  2 


(4.55) 


B=  (5.66  log— 1  +8. 5)2 -2 8. 72  log(h.ie.) 


(4.56) 


Similarly,  substituting  Eq.  4.33  and  4.35  into  Eq.  4.25  and 
4.26,  respectively,  integrating,  and  then  substituting  in 
Eqs.  4.49  and  4.50  provides 


B'~ 


P  si=l  +  - 


6. 
h . 


,,  d-  h,  6.  h  ,,  d„  h ,  6.  _ 

[1O5T*:‘1-04-h;la<>si*:41'04'T)] 


(4.57) 


P 


(4 . 58) 


where 


f?/=  (5.66  log- 


h  Jc. 


-  +8 . 5 ) z-28 


,12  log ( —£h,k,} 
n 


(4.59) 


Graphs  of  0,  /3s1,  and  /3s2  versus  ds/h  are  included  as  Figures 


4.8 


4.12. 
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Figure  4.8  Variation  of  p  for  Logarithmic  Velocity 
Distribution  over  a  Smooth  Boundary 
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Figure  4.9  Variation  of  (3  vs.  ds/h  for  Logarithmic 

Velocity  Distribution  over  a  Smooth  Boundary 
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Figure  4.10  Variation  of  (3  for  Logarithmic  Velocity 
Distribution  over  a  Rough  Boundary 
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Figure  4.11  Variation  of  (3  vs.  ds/h  for  Logarithmic 

Velocity  Distribution  over  a  Rough  Boundary 
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Figure  4.12  Variation  of  /?  vs.  ds/h  for  Logarithmic 

Velocity  Distribution  over  a  Rough  Boundary 
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4.3.3  Power  Velocity  Distribution  for  Turbulent  Flow 


Another  distribution  considered  acceptable  to  represent 
a  turbulent  velocity  distribution  is  a  power  relationship 
such  as 


_i 

m 


(4.60) 


As  before,  the  equation  is  evaluated  and  the  following 
equations  are  revealed 


U-- 


m 

7T7+1 


Uy 


(4.61) 


n  (ml)2 

p  m(m+ 2) 


(4.62) 


P 


Si 


(m+1)2  _n 

m(m+2 ) 


(4.63) 


o  _  (m+ 1)2 
Ps2  m(m+ 2) 


<i  ♦  — ) 
m 


=P( 


u*-) 

m 


(4.64) 


formal,  p  =  =  ps:  =  l .  0208 

d  -2 

Ps?=l  .0208  ( )  6 


(4.65) 


(4.66) 


68 


for  m=l , 


P  =  =  PS_2  =  1 .0159 

d  9/7 

Ps2  =  1.0159  (-£) 


(4.67) 


(4.68) 


A  graph  of  (3,  0s1 ,  and  /?s2  for  m=6  and  m=7  is  given  in  Figure 
4.13. 

As  can  be  seen  in  the  various  graphs,  (3s2  is  not 
particularly  useful  in  this  study.  Originally,  Ps2  was 
developed  in  order  to  possibly  combine  the  effect  of  (3  and 
/?s1 ,  and  therefore  eliminate  one  variable, 
and  )3s1  must  be  used. 


However,  both  (3 
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Distribution 
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5.  DEVELOPMENT  OP  THE  MODEL 

After  identification  of  the  variables  affecting 
sediment  transport,  the  use  of  dimensional  analysis  has 
indicated  which  dimensionless  combinations  of  these 
parameters  are  most  significant.  Theory  has  identified  the 
forces  acting  on  a  single  submerged  particle  and  the  use  of 
the  momentum  equation  has  provided  an  equation  for  sediment 
transport  which  could  be  calibrated  and  verified  given 
sufficient  data.  The  momentum  flux  correction  coefficients 
for  both  the  total  flow  depth  and  that  associated  particle 
depth  have  been  found  in  equation  form  and  are  presented  in 
graphs.  The  employment  of  these  correction  coefficients 
allows  velocity  distribution  to  be  represented  in  the  final 
sediment  transport  equation.  Therefore,  given  the  above 
information,  a  model  is  developed  for  sediment  transport  in 
shallow  flow  on  an  impervious  surface.  Data  is  not 
presently  available  for  this  particular  case  so  assumptions 
are  required  that  when  the  data  is  obtained,  the  chosen 
parameters  can  be  found.  However,  in  order  to  verify  the 
work  of  this  thesis,  some  existing  sediment  transport  work 
is  employed  in  a  modified  manner. 

The  sediment  transport  model  is  developed  by  defining 
an  incipient  motion  parameter  to  predict  particle  movement 
and  an  efficiency  parameter  to  predict  the  portion  of 
particles  set  into  motion  which  actually  are  transported 
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from  the  reach  of  the  channel  under  consideration. 


5.1  INCIPIENT  MOTION  PARAMETER 


As  described  previously,  incipient  motion  occurs  when 
the  forces  acting  to  move  the  submerged  particle  are  equal 
to  those  preventing  movement.  Because  of  the  complexity  and 
nature  of  the  movement  of  a  particle  within  the  flow, 
incipient  motion  is  often  used  as  an  indicator  of  the 
transport.  The  complexities  involved  in  the  study  of 
sediment  transport  prevents  the  determination  of  the 
parameter  from  theory  by  using  the  equations  of  motion. 

Previous  studies  on  the  terminal  fall  velocity  of 
particles  of  different  sizes  often  provides  useful 
information  in  the  study  of  sediment  transport  where  the 
flow  depth  is  much  greater  then  the  particle  size.  However, 
because  of  the  shallow  flow,  the  particle  is  always 
accelerating  or  decelerating  and  never  reaches  its  true 
terminal  velocity. 

However,  it  is  known  that  the  pickup  parameter  is  a 
function  of  the  particle  size  to  depth  ratio,  ds/h,  the 
Reynolds  number  of  the  flow,  Re,  and  the  level  of  turbulence 
in  the  flow.  In  this  study,  the  pickup  parameter  is  defined 
as  when  the  particles  are  set  into  motion  and  is  denoted  by 


n.  n  has  a  value  between  0  and  1. 
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Since  not  all  particles  which  are  set  into  motion  are 
actually  transported  out  of  the  reach  under  consideration,  a 
second  parameter  is  required  for  the  model.  This  parameter 
is  the  deposition  parameter,  M,  and  will  be  related  to  n  by 
an  efficiency  parameter  developed  in  Section  5.2. 

The  pickup  parameter  is  approximated  by  using  the 
Shields  diagram.  Also,  a  second  means  of  finding  the 
parameter  from  known  graphs  of  the  drag  coefficient,  CD, 
versus  Reynolds  number  for  a  sphere  is  discussed. 

5.1.1  Approximation  using  Shields  Diagram 

Figure  5.1  provides  a  modified  Shields  diagram  by  Yalin 
and  Karahan  (1979)  which  can  be  expressed  as 

YCI  =  F(Xcr)  (5.1) 


where 


(5.2) 


Y  =  - - - 

cr  <YS-Y><*, 


(5.2) 


In  addition  to  the  curve  representing  incipient  motion  in 
turbulent  flow,  Yalin  and  Karahan  have  provided  a  curve  for 
laminar  conditions.  The  laminar  line  is  based  on  laboratory 
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(from  Yalin  and  Karahan,  1979) 
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flume  data.  Therefore,  it  is  considered  that  the  two  curves 
can  adequately  describe  both  turbulent  and  laminar  steady 
uniform  flow  which  is  studied  in  this  research. 

The  procedure  for  predicting  motion  given  the  curves  of 
the  modified  Shields  diagram  is  to  find  the  value  of  the 
dimensionless  parameter  Ucr/Jq3s  as  a  function  of  the  Shields 
diagram  and  ds/h,  the  bed  slope  (S0)  ,  and  the  Reynolds 
number  at  which  the  flow  is  capable  of  inducing  particle 
motion,  Reh. 


5. 1.1.1  Steady  Uniform  Flow 


For  a  steady  uniform  flow  r0=YhS0,  the  Shields  diagram 
ordinate  becomes 


(ys-y)de  (ys-y)ds  p -p  d 


(5.4: 


Also  for  steady  uniform  flow,  the  abscissa  is  given  as 


P  -  dsy/?hso  _  y[9at 


TJZ  N 


-fiS0Rec 


(5.5) 


where  Rec  =  Ucr  h/v  is  the  flow  Reynolds  number  at  the 

sediment  incipient  motion  condition. 

The  Shields  diagram  (Figure  5.1)  provides  Eq.  5.4  as  a 
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function  of  Eq.  5.5.  i.e. 


N 


-jf S0Rec=F-' 


[( 


PS“P  d. 


)4-50] 


(5.6) 


Therefore,  the  dimensionless  critical  velocity  relationship 
can  be  expressed  as 


-7=-*\  Ft - i - -£s0] 

N  ^  **  (  P.).! 


(5.7) 


or 


U, 


cr  _ 


’■Re^o  Ft- 


<-^)  -1 
P 


] 


(5.8) 


For  ps/p  =  2.65 


Ucr  _  ,  ds 


■(T)Re^ 


(5.9) 


also  Ucr  -  '  ds 


ds'1  Sc 


^  ^  F[T  Ti?] 


(5.10) 


The  difference  of  the  function  F  for  laminar  and  turbulent 
flows  given  in  Yalin  and  Karahan's  diagram  (Figure  5.1)  is 
due  to  the  difference  in  velocity  distribution  as  reflected 
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by  the  momentum  correction  coefficient  0  as  indicated  by  Eq. 
3.4  from  dimensional  analysis.  Plots  of  Eqs.  5.9  and  5.10 
on  the  basis  of  Yalin  and  Karahan's  diagram  are  provided  as 
Figures  5.2  and  5.3.  Note  that  the  particular  curve 
selected  from  the  Shields  diagram  is  determined  by  the  flow 
Reynolds  number.  Chow  (1959)  suggests  a  transitional 
Reynolds  number  for  open  channel  flow  of  500  to  2,000. 


5. 1.1.2  Steady  Nonuniform  Flow 

As  an  approximation,  the  Shields  diagram  may  also  be 
applied  to  incipient  motion  in  shallow  flow  in  a  wide 
channel  assuming  the  flow  is  under  hydrostatic  pressure 
distribution  (Figure  4.4).  For  a  unit  width,  the  pressure 
difference  between  one  section  (1)  of  the  fluid  volume  to 
another  (2)  can  be  expressed  as 

Ap=-|  y  (hi  ~h\ )  cos0=ycos0  HAh  (5.11) 


The  x  component  of  the  weight  of  the  fluid,  Fhx,  and  the 
shear  force  on  the  boundary,  FT,  of  the  control  volume  are 
expressed  as 


(5.12) 
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Figure  5.2  Critical  Velocity  vs.  ds/h  for  Uniform  Flow 
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Figure  5.3  Critical  Velocity  vs.  ds/h  for  Uniform  Flow 
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Fx  =t0Ax 


(5.13) 


With  incipient  motion  occurring  at  ZFX  =  Fux-FT-Ap  =  0,  the 
following  relationship  is  obtained 

T0=yh  sin0-y  cos0  h  ~  (5.14) 

dx 

Defining  Sf  =  r0/(Yh)  and  noting  that  SQ  =  sin  ©,  we  have 

Sf=SQ-cos0-^  (5.15) 


Thus,  if  the  effect  of  dh/dx=SH  directly  on  the  sphere  is 
assumed  negligible,  the  steady  uniform  relationships  that 
were  obtained  in  Section  5. 1.1.1  can  be  used  as  an 
approximation  provided  S0  is  replaced  by  SQ  -  cos©  SH.  Thus, 
for  small  ©,  from  Eq.  5.9,  noting  that  u.  =  7ghS0, 


^CT  1  -  [  .  ds  S0  Sw 

u.  Rec  h  h  1.65 


(5.16) 


A  graph  of  this  relationship  for  nonuniform  flow  i~  giver,  as 
Figure  5.4. 


5.1.2  Approximation  Using  Particle  Drag  and  Lift  Graphs 


80 


Figure  5.4  Critical  Velocity  vs.  ds/h  for  Nonuniform  Flow 
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Incipient  motion  can  be  induced  by  any  combination  of 
lift,  drag,  or  rotation  as  discussed  in  section  4.1.  The 
case  of  incipient  motion  can  be  expressed  as 

to  [Ft'-F^y)  +(l-o>)  {Fd-Fw)  =0*  (5.17) 

where  <■>  =  a  weighing  factor  for  the  relative  importance 

between  lift  and  drag  on  incipient  motion, 

Fl  is  the  lift  force, 

Fd  is  the  drag  force, 

Fux  is  the  submerged  particle  weight  in  the 
x-direction, 

and  F^  is  the  submerged  particle  weight  in  the  y- 

direction. 

Inserting  Eqs.  4.1,  4.2,  and  4.3  for  the  drag,  lift,  and 
weight,  respectively,  into  Eq.  5.17  yields 

to  [  CL*.U*-±  <  ~  'I )  9ds)  +  ( 1 -to )  [  C^SUZ2 -±  ( -1 )  gdsS0]  =0  (5.18: 
where  SQ  =  sin  ©, 

Ur  =  reference  velocity  at  the  sphere  center, 
cos  9=1  for  a  small  9. 

Since  it  has  been  shown  that  Ucr/Ur  is  a  function  of  (3  and 
/3S,  Eq.  5.16  can  be  rewritten  as 

F(  P  ,  Ps)  [wC^+d-w)  CD]  =4  (  —  -1)  [«+(l-w)50] 
gd '  3  3  p 


(5.19) 
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where  F  =  a  function. 

Values  of  the  lift  and  drag  coefficients  can  be  approximated 
for  a  sphere  from  previously  published  graphs  of  CL  or  CD 
versus  the  flow  Reynolds  number  like  the  one  given  as  Figure 
5.5.  Values  of  a  and  F(/?,/3s)  are  currently  unknown. 

However,  these  parameters  could  be  found  with  sufficient 
experimental  data  when  available.  Once  <o  and  F(/3,0s)  are 
known,  Ucr2/gds  can  be  calculated  since  all  other  variables 
are  known  or  available.  The  critical  velocity  can  then  be 
found  and  compared  agair.st  the  actual  mean  velocity  of  the 
flow  to  see  if  the  flow  is  capable  of  initiating  sediment 
into  motion. 


5.2  SEDIMENT  TRANSPORT  EFFICIENCY  PARAMETER 

It  is  assumed  that  once  incipient  motion  occurs,  the 
flow  has  a  unique  and  specific  transport  efficiency  which  is 
a  function  of  the  local  conditions.  The  sediment  transport 
efficiency  parameter  is  denoted  by  and  it  is  defined  as 
the  fraction  of  the  particles  set  into  motion  which  are 
actually  transported  from  the  reach  under  consideration. 
Since  the  pickup  (n)  and  deposition  (M)  parameters  are 
functions  of  the  same  variables,  the  two  can  be  related  by  a 
new  function, 
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M=K6Q-^  (5.20) 

when  fi  is  greater  than  M,  erosion  occurs  and  when  the 
reverse  is  true,  sedimentation  occurs.  If  we  let  £'=(n-M) 
and  substitute  this  into  Eq.  4.18, 

nV_(Q- Af)=-~£  (5.21) 

s  dx 


where  n  =  the  number  of  particles  available  for 

transport  within  the  given  reach  of  unit  area. 

Therefore,  from  Eq.  5.22 


^-nV,(Q-A;Q-^).nV,Q(l-Xi-^) 


(5.22) 


Therefore,  £ '  =  (n-M)  =  K^,  ranges  from  0  to  1  and  would 
be  determined  through  experimentation,  n,  M,  and  £ '  all 
have  the  units  of  1  per  time  (1//T) .  The  characteristics  of 
can  be  determined  by  examining  fl  and  M.  n  is  a  function 
of  Ucr  as  shown  by 


Qh  _  c,  ,  Ucr 

"5"  'W. 


u 


(5.24) 


Therefore,  n,  M,  and  can  all  be  expressed  as  a  function 
of  the  critical  velocity  expressed  as  (U-Ucr)/Ucr.  Values 
for  the  critical  velocity  expression  is  obtained  from  both 
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physical  theory  and  past  laboratory  work.  The  physical  view 
of  drag,  lift  and  rotation  on  a  submerged  particle  provided 
Eqs.  4.8,  4.11,  and  4.13  where  Ucr  is  expressed  as  a 
function  of  a  K  function  and  d$.  Eqs.  4.8,  4.11,  and  4.13 
can  be  modified  by  the  gravitational  acceleration  constant, 
g,  to  provide 

Drag:  -^SL.  =g  KD  (5.26) 

yfz&s 

Lift :  -^£1—  =g  k  (5.27) 

Rotation:  — — — ■  =g  KD  (5.28) 

where  K  =  a  function. 

The  use  of  Shields  diagram  provided  Eq.  5.10  to  find  a  value 
for  the  critical  velocity  expression. 

A  suggested  plot  of  n  and  M  versus  the  dimensionless 
critical  velocity  expression  (U-Ucr)/Ucr  is  given  in  Figure 
5.6  and  5.7,  respectively.  Both  the  ideal  and  actual  cases 
are  given.  The  actual  case  differs  from  the  real  due  to  the 
bursting  effect  in  turbulent  flow  and  the  temporal 
variations  of  the  wake  in  laminar  flow.  Since  the  velocity 
fluctuates  or  varies  with  time,  the  particle  may  be  set  into 
motion  at  a  particular  instant  when  the  flow  conditions 
warrant  it  even  though  the  average  flow  has  not  reached  the 


Figure  5.6  Pickup  Parameter  vs.  Dimensionless  Critical  Velocity 
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Deposition  Parameter  vs.  Dimensionless  Criti 
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critical  value.  The  deposition  curves  are  negatively  skewed 
because  the  critical  velocity  required  for  pickup  is  greater 
than  that  required  for  deposition. 

Combining  Figure  5.6  and  5.7  provides  a  view  of  the 
efficiency  parameter  since  £'=(n-M).  Figure  5.8  provides 
the  shape  and  characteristics  of  the  efficiency  parameter 
curve.  Since  no  upstream  sediment  concentration  is 
considered,  £'  ranges  from  zero  to  one.  The  must  be 
calibrated  and  verified  by  field  and  laboratory  data  before 
it  can  be  used  practically. 

Since  the  data  is  not  presently  available  to  determine 
1 1 ,  it  is  assumed  that  £  '  is  a  meaningful  parameter  and 
could  be  determined  as  a  function  of  the  variables  affecting 
the  pickup  and  deposition  parameters.  It  is  important  to 
note  that  as  the  flow  becomes  saturated  with  sediment, 
approaches  zero.  Once  the  flow  can  no  longer  support 
additional  sediment  transport,  £'  is  zero  within  that 
saturated  reach.  A  graphical  representation  is  given  in 
Figure  5.9. 

The  procedure  developed  in  this  thesis  is  to  be  applied 
assuming  the  sediment  transport  efficiency  parameter  is 
known  or  can  be  computed.  The  actual  method  used  to  obtain 
£  is  left  to  others  to  develop  in  future  research  from 
either  experimentation  or  theory.  The  theoretical  procedure 
would  be  to  obtain  £  directly  from  the  n  -  M  relationship 
since  both  all  are  a  function  of  the  same  variables.  In  a 
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Figure  5.9  Variation  of  £'  with 
Sediment  Concentration 


0  i - 

0  C,  concentration 


laboratory  setting,  £  can  be  obtained  by  using  flume  data 
where  £  is  obtained  as  a  function  of  time  or  space  depending 
upon  the  arrangement  of  the  flume  and  the  data  recorded. 

Once  the  £  relationship  is  developed,  the  transport 
efficiency  parameter  can  be  applied  to  different  situations 
and  conditions. 

For  ease  of  use,  a  dimensionless  efficiency  parameter 
is  given  by 


c_  (Q  -AD 

*  Q 


(5.25) 


5.3  SEDIMENT  TRANSPORT  EQUATION 

The  transport  rate  can  be  approximated  by  using  an 
assumed  transport  efficiency,  £ ,  and  given  a  procedure  to 
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determine  whether  the  flow  is  sufficient  to  induce  particle 
motion.  The  sediment  transport  term,  qs,  is  defined  as  the 
total  sediment  volume  transported  per  unit  width  of  the 
channel  and  per  unit  time.  From  the  continuity  equation 
shown  previously,  the  general  form  of  the  equation  is 

+  (5.25) 


where  n  =  number  of  particles  per  unit  bed  area, 

i  =  the  transport  efficiency  parameter, 

Vs  =  volume  of  a  sediment  particle, 

Ax  =  length  of  the  channel  segment. 

A  plot  of  Ts/q  td  vs.  time/td  is  given  as  Figure  5.10. 
Since  the  Eq.  5.23  represents  the  sediment  transport  rate 
for  a  unit  time,  the  total  transport  for  an  entire  event  of 
duration  td,  Ts,  over  a  length  x=0  to  x=Lg  with  qsQ=  0,  is 

td  A) 

Tg=f(fV(x,t)  nix.t)  Vsdx)dt  (5.29) 

0  0 

For  Lq  =  E  Ax  and  duration  td  =  Z  At  with  £  '  invariant  with 
time  but  variable  with  space, 

r.-EE 

i-i  j-l 


(5.29) 
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Figure  5.10  Total  Sediment  Transport  vs.  Time 
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since  for  any  given  Ax  the  amount  of  sediment  on  the  bed  at 
the  ith  At  is  (l-£')nVs  and  the  sediment  transport  is  £'(l- 
£')i'1nVs.  Figure  5.11  provides  a  plot  of  T8lt/Tsltd  vs. 
time/td  where  Tslt  is  the  total  sediment  transport  per  unit 
time  and  Tsltd  is  the  total  sediment  transport  per  the  entire 
event.  Eq.  5.25  expresses  the  total  sediment  transport,  T  , 
as  a  volume  per  unit  width.  Although  mass  is  easier  to 
measure  and  more  commonly  used,  volume  provides  a  suitable 
variable  from  the  physical  viewpoint.  Also,  mass  can  be 
easily  found  by  knowing  Vs  and  ps  or  ys*  As  mentioned 
previously  in  Section  5.2  and  Figure  5.8,  £’  decreases  as 
the  sediment  concentration  increases.  Therefore,  Eq.  5.25 
(finite  one)  is  valid  only  for  low  concentrations  where 
is  approximately  constant.  Figure  5.12  shows  a  graph  of  the 
dimensionless  parameters  qs/q  versus  x/h  computed  using 


<?sx_  Vnvs  x 

g  U  h 


(5.30) 


However,  the  transport  efficiency  parameter  is  most 
likely  not  constant  throughout  the  channel  reach.  As  the 
flow  becomes  saturated  with  sediment  and  is  no  longer  able 
to  induce  additional  transport,  the  efficiency  should 
decrease  to  zero.  The  decrease  of  £  must  be  determined  by 
data  analysis.  However,  an  assumed  relationship  is  used  in 
this  study.  Figure  5.13  graphs  qs/q  versus  x/Lg  for  a 
decaying  i 1 .  The  decay  relationship  is  assumed  to  be  an 
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Figure  5.11  Total  Sediment  Transport  vs.  Time 
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Figure  5.12  Plot  of  qs/q  vs.  x/h 
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Figure  5.13  Plot  of  qs/q  vs.  x/L0  with  Decreasing 
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exponential  one  of  the  form 


e-:<x 


(5.31) 


where  k  is  the  decay  rate  (1/length). 


Figure  5.14  qs/qs  mai(  versus  Sediment 
Concentration 


The  sediment  transport  capacity  of  the  flow  may  never 
approach  saturation  if:  the  efficiency  is  constant;  the 
sediment  concentration  is  very  sparse;  or  the  channel  length 
is  comparatively  short.  See  Figure  5.14.  Under  one  or  more 
of  the  above  conditions,  the  plot  of  gs/q  under 
consideration  will  be  a  straight  line  instead  of  the 
presumed  line  asymptotic  to  the  saturation  level.  Such  may 
be  the  case  on  an  airfield  pavement  where  the  channel  is 
relatively  short  (i.e.  across  the  runway  or  apron). 
Therefore,  the  estimation  of  a  constant  £  may  be  valid  for 
this  particular  case. 
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6.  RESULTS  &  DISCUSSION 


6 . 1  RESULTS 


Since  field  or  laboratory  data  is  not  available  for 
sediment  transport  in  shallow  flow,  hypothetical  data  is 
used  to  test  the  model  developed  in  the  previous  section. 
Figure  5.12  shows  a  plot  of  qs/q  versus  x/h.  The  shape  of 
the  curve  and  the  magnitude  of  the  values  seem  reasonable 
given  the  physical  view  of  the  process  and  the  selected 
values  of  the  variables. 

Although  a  lack  of  available  data  prevents  a  calibrated 
and  verified  quantitative  model,  a  qualitative  one  is 
provided  by  the  research  in  this  thesis.  This  model  can  be 
applied  to  simulate  sediment  transport  in  shallow  steady 
uniform  flow  over  an  impervious  surface.  Future  expansion 
of  the  results  presented  here  will  allow  the  model  to  be 
used  for  nonuniform  and  unsteady  flows  also.  Therefore,  the 
objective  of  the  thesis  has  been  obtained  in  the  form  of  a 
model  which  predicts  the  total  sediment  transport  in  a 
shallow  flow  over  an  impervious  surface. 
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6.2  DISCUSSION 


6.2.1  Sediment  Transport  Parameters 

The  major  parameters  presented  in  this  study  are  the 
incipient  motion  parameter,  Ucr/./gc[s,  and  the  transport 
efficiency  parameter,  £'.  The  use  of  dimensional  analysis 
and  the  approach  from  general  physics  helps  to  verify  the 
validity  of  these  two  parameters.  The  employment  of  the 
Boussinesq  type  coefficients  ensure  the  incipient  motion 
parameter  accounts  for  the  particular  velocity  distribution 
(laminar  or  turbulent)  and  the  small  flow  depth  to  sediment 
size  ratio.  Unfortunately,  a  value  for  the  transport 
efficiency  parameter  cannot  be  obtained  due  to  the  lack  of 
data.  However,  this  parameter  is  valid  and  the  model  can 
account  for  variation  of  $ .  An  expression  for  $  as  a 
function  of  the  other  sediment  transport  variables  was  not 
developed  due  to  the  complex  nature  of  the  flow  and  the  lack 
of  data.  Although  the  efficiency  is  assumed  constant  over 
the  entire  length  of  the  channel  under  consideration,  a 
constant  £  may  be  applicable  to  the  actual  case  of  sediment 
transport  in  shallow  flow  over  an  impervious  surface  such  as 


an  airfield. 
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6.2.2  Impact  of  Variables  Not  Considered 

Due  to  the  complexity  of  the  problem  of  particle 
movement  in  a  moving  fluid,  many  variables  have  not  been 
considered  as  a  separate  parameter  in  this  study.  Unsteady 
flow,  the  impact  of  rainfall,  and  the  effect  of  mutual 
interference  of  multiple  particles  are  three  such  additional 
parameters.  The  impact  and  approach  to  inclusion  of  each 
parameter  is  discussed  below. 

Unsteady  flow  can  be  included  in  the  sediment  transport 
parameter  by  modeling  the  velocity  distribution.  The  same 
general  form  of  the  equations  presented  earlier  are  still 
valid.  However,  the  coefficients,  such  as  /3,  will  be 
different.  Also,  in  order  to  account  for  unsteady  flow,  the 
model  would  require  a  great  deal  more  data.  One  example  is 
the  variation  of  the  velocity  with  time  which  would  be 
required  to  replace  the  average  velocity.  This  requirement 
for  additional  data  would  most  likely  not  be  greater  than 
the  benefit  of  a  more  accurate  model. 

The  effect  of  rainfall  on  sediment  transport  on 
erodible  boundaries  has  been  studied  in  recent  years. 
Generally,  it  is  presented  that  the  impact  of  the  raindrops 
imparts  energy  to  the  control  volume.  The  friction  slope 
term,  Sf,  is  often  modified  to  include  the  effects  of  the 
raindrop  impact.  However,  no  work  has  been  found  which 
studied  rainfall  impact  on  sediment  transport  over  an 
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impervious,  nonerodible  boundary.  Much  laboratory  flume 
work  is  required  to  establish  the  necessary  relationships  of 
the  raindrop  variables  with  the  sediment  transport  variables 
presented  here  earlier. 

Only  a  single  particle  was  considered  in  the  earlier 
theory  section.  Figure  6.1  depicts  the  forces  acting  on 
submerged,  adjacent  particles.  In  order  to  relate  these 
forces,  a  new  term  is  introduced.  The  term  Fp  denotes  the 
"pushing"  force  given  to  one  particle  by  another.  As  seen 
in  the  figure,  the  major  impedance  to  the  movement  of 
particle  number  2  is  the  shielding  it  receives  from  particle 
1.  This  shielding  prevents  particle  2  from  receiving  the 
full  impact  of  the  flow.  As  stated  before,  the  term  Asn, 
particle  surface  area  normal  to  flow,  can  partially  account 
for  the  shielding.  However,  this  term  does  not  account  for 
the  variation  of  flow  caused  by  particle  1,  i.e.  the  flow 
separation  and  resulting  low  pressure  area  which  will 
develop  immediately  behind  particle  1.  This  pressure 
reduction  results  in  greater  drag  on  particle  1  but  will 
tend  to  decrease  the  total  drag  on  particle  2. 

This  decrease  in  total  drag  on  particle  2  is  due  to  the 
decreased  (actually,  negative  according  to  the  defined 
coordinates)  pressure  drag  on  the  particle.  The  pressure 
drag  and  skin  friction  combine  to  account  for  the  total 
drag.  The  term  C0,  drag  coefficient,  accounts  for  both 
types.  This  effect  of  reduced  drag  on  the  shielded  particle 
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will  vary  greatly  according  to  the  flow  regime,  laminar  to 
fully  developed  turbulent.  At  low  Reynolds  number,  Re  <  0.1 
-  1.0,  the  flow  is  considered  creeping  and  skin  friction 
drag  is  predominate  (Shames,  1982) .  Therefore,  the  effect 
is  minimal.  However,  as  Re  increases.  Re  ~  5xl03, 
separation  occurs  and  pressure  drag  becomes  important 
(Shames,  1982) .  When  the  flow  becomes  increasingly 
turbulent,  the  pressure  drag  is  reduced  because  the 
turbulent  flow  is  better  able  to  adhere  to  the  curved 
surface  of  the  particle  and  resist  separation. 

When  the  pressure  behind  particle  1  is  reduced, 
particle  2  will  tend  to  move  toward  particle  1.  This  is  due 
to  the  larger  pressure  behind  particle  2  (supposing  no  other 
particles  are  close  to  particle  2) .  However,  the  movement 
of  the  two  particles  toward  each  other  will  increase  the 
fluid  pressure  between  the  two  particles.  When  this 
increased  pressure  becomes  greater  than  the  pressure  behind 
particle  2,  the  particles  will  begin  to  move  away  from  one 
another.  This  action  of  particle  movement  toward  and  away 
from  each  other  can  be  accounted  for  by  the  "pushing"  force 
term,  Fp.  However,  details  one  the  magnitude  or  makeup  of 
this  term  are  not  yet  known. 

A  conspicuous  assumption  used  in  the  model  presented 
here  is  that  the  sediment  transport  rate  can  be  found  by 
relating  it  to  conditions  of  incipient  motion  rather  than 
modeling  the  sustained  movement  of  the  particle.  However, 
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this  approach  is  the  one  generally  taken  by  those  studying 
sediment  transport  in  alluvial  cases.  The  complexities  and 
lack  of  greater  physical  understanding  of  the  equation  of 
motion  for  the  particles  virtually  require  such  an  approach. 

6.2.3  Application  of  the  Model 

The  sediment  transport  procedure  developed  in  this 
thesis  can  be  applied  given  values  for  the  required  data. 
Sediment  variables  include  ds,  ps,  nVs,  and  £.  The  first  two 
are  found  from  analysis  of  the  field  data  sample.  The  term 
nVs  can  be  used  to  replace  Ap,  the  areal  sediment  packing, 
and  expresses  the  total  volume  of  sediment  available  for 
transport.  Since  the  particle  diameter  is  known,  the 
particle  volume  is  easily  calculated  for  the  case  of  a 
sphere.  The  volume  of  a  non-spherical  particle  could  be 
related  to  that  of  a  spherical  one  by  the  use  of  the  shape 
factor,  SF.  Since  the  total  volume  is  sought,  the  actual 
number  of  particles,  n,  is  not  required.  This  total  volume 
variable,  nVs,  could  be  obtained  by  statistical  methods 
which  relate  the  total  available  volume  to  factors  such  as 
location,  rain-free  days,  frequency  of  street  cleaning,  etc. 
for  urban  environments.  It  must  be  noted  however,  that  the 
airport  environment  greatly  differs  from  that  of  the  urban 
one,  mainly  because  the  jet  engine  blast  from  the  aircraft 
will  continuously  move  the  sediment.  The  actual  method 
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employed  to  determine  the  available  source  of  sediment  is  a 
problem  not  addressed  in  this  work. 

Fluid  variables  include  p,  q,  h,  /i,  y,  and  For 

rainfall,  p,  fi,  and  y  are  fairly  constant  and  known.  The 
flow  rate  q,  and  flow  depth  h  can  be  measured  directly  in 
the  field.  Given  the  geometry  of  the  flow,  either  q  or  h 
could  be  replaced  by  U,  the  mean  flow  velocity  since  q  =  hU. 
/3  is  provided  in  Section  4.3.  F.1  ow  variables  are  S0  and  L0 

which  can  be  found  directly  from  the  field. 

Therefore,  the  only  variable  above  which  is  not  easily 
found  or  expected  to  be  found  with  reasonable  means  is  £ , 
the  sediment  transport  efficiency  parameter.  The  procedure 
to  determine  this  variable  is  left  to  others  as  outlined  in 
Section  5.3. 

Particles  of  different  sizes  or  densities  can  be 
accounted  for  in  the  model  by  classifying  the  different 
sizes  or  densities  into  groups,  computing  the  transport  for 
each  group  and  then  adding  the  transport  values  together  for 
a  total.  This  assumes  a  linear  approximation  is 
appropriate.  Further  study  is  required  to  verify  this 
assumption. 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  physically  based  procedure  to  simulate  sediment 
transport  is  presented  for  the  case  of  shallow  flow  over  an 
impervious  surface.  The  procedure  is  developed  for 
particulate  matter  only  and  would  not  apply  to  oils,  fibrous 
material,  or  dissolved  pollutants.  This  model  accounts  for 
either  a  laminar  flow  following  a  parabolic  distribution  or 
turbulent  flow  over  a  smooth  or  rough  boundary  following  a 
logarithmic  distribution  or  a  power  distribution.  Momentum 
flux  correction  factors  and  established  are  presented  to 
account  for  the  variable  velocity  distribution  and  varying 
flow  depth  to  particle  size  ratio. 

No  field  or  laboratory  flume  data  was  collected  as  part 
of  this  research.  Such  data  was  not  found  available  in  past 
research  by  others  although  this  particular  case  of  sediment 
transport  seems  to  be  gaining  interest  in  the  academic 
environment.  Although  the  model  is  not  verified  by  actual 
data,  hypothetical  data  is  used  to  show  the  relationships 
developed  herein.  All  graphs  contained  the  shapes  and 
characteristics  expected  from  physical  knowledge. 

Because  of  this  lack  of  calibration  and  verification, 
the  model  should  not  be  applied  to  any  actual  case  for 
prediction  purposes  until  tested.  As  Bagnold  (1966)  stated: 
"Indeed  the  outlook  is  obscured  by  the  very  number 
of  published  theories.  For  example,  the  engineer 
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research  student,  encouraged  by  degree  thesis 
requirements  to  advance  a  new  theory  rather  than 
to  verify  existing  concepts,  is  able  to  find 
somewhere  in  the  literature,  and  to  quote,  support 
for  an  assumption  which  best  suites  his  argument." 

Usually,  a  sensitivity  analysis  is  conducted  to  reflect 
the  sensitivity  of  the  model  to  each  parameter.  Therefore, 
the  parameters,  when  erroneous,  are  identified  which  impart 
the  most  error  in  the  model.  No  sensitivity  analysis  was 
accomplished  in  this  thesis  since  the  data  is  not  available 
to  fully  test  the  model.  Since  the  data  must  have 
sufficient  accuracy  in  order  to  have  a  meaningful 
uncertainty  range,  a  sensitivity  analysis  would  not  be 
helpful  in  this  work.  However,  the  graphs  presented  in 
Section  5.3  are  given  as  functions  of  dimensionless 
parameters  in  order  to  reduce  any  sensitivity  of  the  model 
to  a  particular  variable. 

The  major  recommendation  for  future  study  is  the 
collection  of  applicable  laboratory  flume  and  field  data. 
This  data  is  necessary  to  verify  the  relationships 
established  in  this  study  and  provide  means  to  expand  and 
detail  these  relationships.  Section  6.2.2  discusses  several 
variables  which  were  not  fully  incorporated  into  the 
procedure  developed  in  this  thesis.  The  expansion  of  the 
model  to  fully  incorporate  these  variables  is  a  requirement 
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of  future  study.  These  variables  include:  rainfall; 
particles  of  varying  shapes,  sizes,  and  densities;  and 
nonuniform  and  unsteady  flow.  Future  study  is  required  to 
develop  the  sediment  transport  efficiency  parameter.  Also, 
a  method  to  determine  the  total  available  sediment  is 
needed. 
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